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1 Introduction, General Idea and Outline 
1.1 Introduction and outline 

The AdS/CFT conjecture originally proposed by Maldacena [1], refined in [21 [3] and reviewed 
in [1] has been one of the most interesting developments in theoretical physics of the last 
decades. It has become one of the most powerful analytic tools to deal with strong coupling 
effects of some particular gauge theories in the planar limit A^ — ?■ 00. The most studied and 
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best understood case corresponds to SU{Nc) A/" = 4 SYM which is a highly supersymmetric 
conformal theory and which only contains matter in the adjoint representation of the gauge 
group. Certainly, there are many interesting field theories which do not share these properties 
and this fact has lead to an enormous amount of effort devoted to extending the duality along 
different paths. Consequently, people have constructed gravity duals of non-supersymmetric, 
non-conformal gauge theories, in different vacua and with diverse matter contents. One 
can mention the attempt of building a dual as close as possible to QCD as an aim for 
these generalizations. However, one should keep in mind that this is just one among many 
desirable motivations, since understanding gauge theories at strong coupling (or using gauge 
theories to understand gravity) is a very relevant problem per se, with both theoretical and 
phenomenological possible implications. 

An important development of AdS / CFT has been to generalise the matter content of the 
gauge theories under consideration and, in particular, to include fields which transform in the 
fundamental representation of the gauge group, as the QCD quarks dc0. A first possibility is 
to add the fiavors in the quenched approximation. The word quenched comes from the lattice 
literature and, in that context, it amounts to setting the quark fermion determinant to one. 
In more physical terms, quenching corresponds to discarding the quark dynamical effects. 
This means that quantum effects produced by the fundamentals are neglected, the quarks 
are considered as external non-dynamical objects in the sense that they do not run in the 



loopqj. From the string side, adding quenched quarks to a given gauge theory corresponds 
to incorporating a set of brane probes in the dual background, which is not modified with 
respect to the quark-less case. By analysing the worldvolume physics of these fiavor branes 
(typically using the Dirac-Born-Infeld -|- Wess-Zumino action) a lot of physically interesting 
questions can be understood. For instance and just to mention a few, chiral symmetry 
breaking can be neatly described, phase diagrams can be constructed and meson spectra can 
be exactly computed. It is hard to do justice to the huge literature in the subject so let us 
just mention the seminal papers |5], [6] and a recent review [7]. 

Thus, it is fair to say that the study of quenched fiavor within the gauge-gravity corre- 
spondence has been a very fruitful program. Nonetheless, there are physical features which 
are intimately related to the quantum effects of the quarks. Examples are the consequences 
of the presence of fundamentals on the running couplings, which may ultimately lead to 
conformal points, conformal windows [8] or Seiberg-like dualities [9]. More phenomenologi- 
cally, multihadron production, the screening of the color charge or the large mass of the t]' 
meson are spin-offs of these quantum effects. Let us also mention that the most successful 

"'^With an abuse of language, we will use throughout this review the words quark or flavor to refer to any 
field, fermionic or bosonic, transforming in the fundamental representation of the gauge group. Accordingly, 
by mesons we will mean bound states of quarks. 

^In the lattice, usually, quenching is thought to be a good approximation for heavy quarks whereas for 
the gauge-gravity examples the relevance of the quenched approximation comes from having parametrically 
less fundamental than adjoint fields Nf <C Nc. 
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application of string duals towards phenomenology has been the construction of solutions 
that can be used as toy models for the experimental quark-gluon plasma. Thus, a very 
interesting program is to build black hole solutions with unquenched flavor which are really 
dual to quark-gluon plasmas, i.e. such that the effect of the dynamical quarks affects the 
plasma physics, as is expected to be the case in the real world. 

These observations largely motivate the study of theories with unquenched quarks from 
the string theory dual point of view. Unquenching the flavors of the gauge theory has a 
very precise implication for the dual theory: the gravity background has to be modified by 
the inclusion of the quarks, namely, one needs to take into account the back-reaction on the 
geometry produced by the flavor branes. The main goal in the following will be to present 
methods to compute such back-reacted solutions. This will be done by presenting different 
examples that, hopefully, will help the reader to gain insight in both the physical questions 
and the technical tools used to address them. 

In this review, we will focus on a specific family of unquenched constructions. Namely, 
we will discuss at length just solutions of type IIA or type IIB string theory in which the 
fundamentals come from a smeared set of flavor branes. In section II. 5[ we will try to provide 
a general understanding of this notion of smearing the flavor and argue why we find it a case 
of particular interest. As we will see, by considering the case of smeared D-branes we can 
build a systematic approach applicable to different situations and which typically results in 
large simplifications as compared to other kind of flavor D-brane distributions. This smearing 
procedure referring to flavors was first introduced in [10] in a non-critical string framework 
and in [11] in a well-controlled ten-dimensional context. 

It is important to remark that this smearing is by no means the only possibility to introduce 
unquenched fundamentals in gauge-gravity duals. Many important works have followed 
alternative paths to construct different models. We are not able to review them here, but 
we provide a survey of the literature in section [L6l 

Outline 

We will devote the rest of section [1] to further clarifying the kind of physical problems we 
want to address and to give the general methods and notions which are common to all the 
constructions we will present later. 

Then, sections [2ll6] will analyse different models that are ordered in increasing order of 
complexity. Each section can be read mostly independently from the rest. The discussion of 
each model can always be regarded as a two step process. First, one has to solve the equations 
for finding the back-reacted solutions of type II supergravity coupled to a set of D-brane 
sources. Second, one can use these solutions to extract the physics of the conjectured gauge 
theory duals with unquenched flavors. Readers interested in different aspects of the problem 
can consult the different parts independently. We would like to stress that, even without 
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making any reference to the gauge-gravity correspondence, the string theory solutions and 
methods developed to find them are interesting by themselves. 

Section |2] deals with the backreaction of D7-branes on AdS^ x spaces, where stands 
for a Sasaki-Einstein space. As a matter of fact, a large part of the discussion can be carried 
out without specifying the X^. Notwithstanding, the two most interesting cases correspond 
to X^ = and X^ = T^'^. At different points during section [21 we will refer to these 
particular examples in order to explain concrete features. We will present supersymmetric 
solutions for massless and massive quarks, and also non-supersymmetric black hole solutions 
which are dual to theories at finite temperature, in a deconfined plasma phase. All these 
solutions share the property of having a singularity, associated with a UV Landau pole 
in the field theory (when quarks are massless and the temperature is zero, there is also 
a naked IR singularity). We will show how to make well-defined IR predictions from the 
geometry, even in the presence of the UV singularity (in much the same spirit as in field 
theory renormalization). 

In section [31 we will discuss a model in which both the color and fiavor branes are D5's. 
It is dual to a (3+l)-dimensional A/" = 1 theory with a UV completion. Among several 
nice features of the model that will be presented, we would like to remark here that it 
incorporates a geometrical description of a Seiberg-like duality. Section [H is also built from 
a D5-D5 intersection and in fact shares several similarities with the previous model. The 
construction corresponds to color D5's wrapping a compact 3-cycle and therefore the dual 
field theory is (2-1-1 )-dimensional. 

In section [5] we examine the addition of D7-branes to the conspicuous Klebanov-Strassler 
model [12]. From the physical point of view, how the unquenched fiavors affect a duality 
cascade is particularly interesting. From a technical point of view, the system is slightly 
more involved than the previous ones because different RR and NSNS forms are turned on. 
However, despite this complication, it is remarkable that almost all functions of the ansatz 
can be integrated in a closed form. 

Section [6] reviews a different class of models. The dual gauge theories are built on wrapped 
supersymmetric D-branes with the peculiarity that some of the adjoint scalars remain mass- 
less. As we will explain, it is not sensible to smear the branes in all the transverse directions. 
The associated main technical difficulty will be the fact that one has to solve partial differ- 
ential equations to find the background. 

Profiting from the experience gained by discussing these examples one by one, in section [71 
we will give a more mathematical viewpoint of the constructions. In particular, we will take 
some tools of differential geometry to describe in a concise and compact way the distributions 
of mass and charge due to the presence of the flavor branes. 

Finally, in section [8] we will conclude by summarizing the whole topic trying to give a 
general perspective of the results obtained and by also providing an outlook of the subject. 
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1.2 Presentation of the problem 

As anticipated in the introduction, we will discuss the addition of flavors to field theories 
(mostly focusing on SUSY examples, but this is not mandatory) using AdS/CFT or more 
generally, gauge-strings duality. 

We hope it is clear to the reader that the addition of flavors (fields transforming in the 
fundamental representation) is a very interesting exercise from a dynamical point of view. 
Indeed, in a theory with adjoint fields (let us for the sake of this discussion, consider the case 
of a confining field theory) the presence of fields transforming in the fundamental will produce 
the breaking of the "QCD-string" or screening. Of course, the fundamentals will add a new 
symmetry, that can be SU{Nf) or, like in massless QCD, SU{Nf) x SU{Nf); a baryonic 
U{1)b symmetry should also appear. Obviously the presence of global symmetries (and their 
possible spontaneous or explicit breaking) will directly reflect in the spectrum. Apart from 
this, it will happen that the states before the addition of fundamentals, that is the glueballs, 
will interact and mix with the mesons, giving place to new diagonal combinations that will 
be the observed states. Moreover, anomalies will be modified, as fermions that transform 
in the fundamental will run in the triangles. Also gauge couplings will run differently and 
finally new dualities (Seiberg-like [9]) may appear. In the rest of this article, we will discuss 
how all of the above mentioned features are encoded in string backgrounds. 

It is clear that we need to add new objects to our string background. These new objects 
are D-branes, on which a gauge field propagates, encoding the presence of a U{Nf) gauge 
symmetry (in the bulk), dual to the global U{Nf) in the dual QFT. Also, it is on these D- 
branes that the meson fields, represented by excitations of the branes, propagate and interact. 
Following a nomenclature that by now became standard, we will call these D-branes "flavor 
branes" . 

It is then clear that to add flavors to a field theory whose dual we know, we should consider 
the original (unflavored) string background and add flavor-branes. Now, the point is how to 
proceed technically to add these new branes. It may be useful to consider two developments 
of the 1970's, that will turn to illuminate on the answer to this question. 

In the papers |13j and [H] 't Hooft and Veneziano respectively considered the influence of 
fundamentals when the following scaling is taken 

oo, A = Qym^c = fixed, (1.1) 

and considered the two possible cases ('t Hooft and Veneziano respectively) 

Nf 

Nf = fixed, X = — )■ 0, 
Nf 

A^/ — )■ oo, X = — = fixed . 
It is very illuminating to see how different diagrams contributing to the same physical process 
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(to fix ideas, an n-point correlator of mesonic currents) scale in these two cases. In this 
respect, a formula for the kinematical factor of the scattering of n mesons was produced in 
|15j . considering diagrams with w internal fermion loops (windows), h non-planar handles 
and h boundaries, 

<B,...B„>~(^)'°ivf--^'-". (1.2) 

Consider the case of scattering of two mesons n = 2. We see that diagrams like the first one 
in the figure (w = /i = 0, 6 = 1, n = 2) scale like a constant ~ 1, the second diagram (with 
iy = l,/i = 0, 6=1,^ = 2) scales like while the third one (with = 0, = 0, 6 = 2, n = 2, 
that is non-planar) goes like N~^. 





Figure 1: Diagrams for a meson propagator, with two insertions of the meson operator [n = 2) 
shown as thick points on the boundaries. The dashed lines are gluons that fill the diagram in the 
large Nc limit and the thick lines are quarks. A) Planar diagram with no internal quark loops 
(/i = 0, li; = 0), the scaling is ~ 1. B) Planar diagram with an internal quark loop /i = 0, w = 1, 
the scaling is ~ Nf/Nc- C) Non-planar diagram with no internal quark loops /i = 0, u; = 0, 6 = 2, 
the scaling is ~ l/Nc- 

So, we see that from this view point, the Veneziano scaling captures more physics, rep- 
resented here by diagrams like B in Fig. [H Nevertheless, there may be some particular 
problems for which studying things in the 't Hooft scaling may be enough. 

From the view point of a lattice theorist, working in the 't Hooft scaling, hence neglecting 
the effects of fundamentals running inside loops, is the same as working in what they would 
call the "quenched approximation". We can think of the field theory as being quenched 
when the fundamental fields do not propagate inside the loops. One natural way of doing 
this is to consider the case of very massive quarks. Indeed, when quenching, one considers 
an expansion of the fermionic determinant (for massive quarks) of the form: 

logdet[7''L'^ + m]= Tr\og{j^D^ + m) = Tr log(m[l + ^]) ~ Tr log(m) + 0(l/m) . 
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Keeping only the constant term (or considering a very large mass) is equivalent to saying 
that fundamentals are very difficult to pair-produce, hence their presence inside loops will be 
very suppressed. Another way to quench in the field theory is to consider the case in which 
the quotient x = ^ is very small. Notice that the quenched theory is not equivalent to a 
theory with only adjoints, as fundamentals can occur in external lines, like in a correlator 
of two mesonic currents as exemplified in the diagram A of FigJH Needless to say, lattice 
theorists developed techniques to quench fundamental fields with arbitrary mass. Also, while 
at first sight the quenching as described above is not a good operation as it breaks unitarity 
(not including all possible diagrams) this kind of troubles will be avoided when working in 
the 't Hooft scaling, where unitarity problems will be suppressed in 1/Nc (but of course will 
be present in a lattice version of theories with finite Nc). 

The interesting point to be taken from this, by a physicist working in gauge-string duality 
is that both scalings ('t Hooft's and Veneziano's) can be realized with D-branes. Indeed, in 
both cases we must add D-branes (to realize symmetries and new states as discussed above) , 
but we can add these fiavor branes in two ways: 

• We can add Nj fiavor branes in such a way that we will only prohe the geometry 
produced by the N^. color branes. In this case the dynamics of the probe-fiavor branes 
(the mesons) will be infiuenced by the dynamics of the color branes (the glueballs) 
but not vice-versa. This is a good approximation if x = ^ — )■ 0, which immediately 
sets us in the 't Hooft scaling limit. Notice however, that when the A^^^ contribution 
to some particular quantity vanishes, the fiavor effects may be the leading ones even 
when Nf <t: Nc. 

• We can add Nf flavor branes, in such a way that they will deform the already existing 
geometry, in other words backreacting on the original "color" geometry. In fleld theory 
language, we would say that the dynamics of the glueballs and that of the mesons 
influence each other, leading to new states that will be a mixture of mesons and glue- 
balls. This is surely what we need to do if a; = ^ = fixed and doing this will set us 
in the Veneziano scaling limit. 

More technically, in the 't Hooft scaling limit we are studying the Born-Infeld-Wess-Zumino 
dynamics for a Dk flavor brane in a background created by Nc color Dp-branes (we will 
always work in Einstein frame in the following): 

Sbiwz = -Tk J d'+'xe'^'^^det[gab + 2na'e-lj^ab] + ^ J C A , (1.4) 

where gab, J^ab are flelds induced by the color branes background on the (few) flavor branes. 
The 'shape' of the flavor branes (induced metric) will then influence the mass spectrum 
and interactions of the fluctuations of the flavor branes (the mesons), explicitly realizing 
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the picture advocated above. This hne of research was initiated by Karch and Katz [5] and 
substantially clarified in subsequent papers [6], [16]- [Sn]. This line kept on growing in the 
last few years, finding numerous applications. See the paper for a comprehensive review. 

On the other hand, working in the Veneziano scaling limit implies that we will need to 
study the action 

S = SiiA/IIB + Sbiwz- (1-5) 

There will be new equations of motion, encoding explicitly the numbers Nf, N^. As discussed 
above, it is now very clear that proceeding like this will be the only possibility when the 
number of fiavors is comparable with the number of colors. Also, it makes manifest the fact 
that the dynamics of glueballs (represented on the string side by Sjja/iib which is of order 
g^^ ~ A^"^) is infiuenced and infiuences back on the dynamics of fundamentals (represented 
by Sbiwz, of order gJ^Nf ~ NcNf). The rest of this review will focus on this second scaling 
(Veneziano) . 

Notice that (in both scalings) we are making an explicit difference between the color 
SU{Nc) gauged symmetry and the flavor SU{Nf) global symmetry on the fleld theory side. 
From the string theory construction, this qualitative difference is connected to the fact that 
the volume of the flavor branes is inflnite, as compared to the volume of the color branes. 
Indeed, in the bulk, we only need to realize the fleld theory global symmetry, and we do it 
with the gauge fleld present in the Born-Infeld-Wess-Zumino action. Searching for solutions 
of Dp color branes in interaction with Dk flavor branes in pure IIA/IIB supergravity is an 
interesting problem, but will not represent the physical system we are after, as only flavor 
singlet states would be included in the dynamics. 

Before we proceed studying the formalism and examples to clarify the details, some com- 
ments are in order. 



1.3 The string action and the scaUng hmit in Nc and Nf 

Let us study a bit more the expression of eq. (11. 5p . being careful about coefficients. We will 
consider the case of a set of Nc "color" Dp-branes and Nf "fiavor" Dfc-branes. The action 
for this system will be, in Einstein frame: 



S 



^6 

2 X /! 



S 



-NfTk [ d'+^x e^'t>\l det[gab + 27ra'e-t J",,] + NfT^ [ 

J Jk 

L{IIA/IIB) - 2KloNfTk [ Lbiwz 



+ J CS — terms 



2^10 



;i.6) 
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where by Fi we have denoted the various RR fields and with CS — terms the possible Chern- 
Simons terms. We have taken the simplification of writing the action for the set of fiavor 
branes as Nf times that of a single D-brane, which is enough for the large counting we 
want to undertake here. The gravitational constant and D-brane tension are: 

2k\, = (27^)\2^'^ n = \ , . (1.7) 
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The typical quantization condition for the color branes reads 

1 



2^10 

As a consequence of (11.61) . we will have equations of motion, that generically will read for 
the metric and dilaton: 



L[IIA/B] - 2Kif^NfTkL[BIWZ] , (1.9) 



together with the modified (by the CS-terms) Maxwell equations and importantly, the 
Bianchi identity for the (magnetic) Ramond-Ramond field -Fs-fc that couples to the fiavor 
Dfc-branes: 

dFs-k = 2Kl^NfTk5''-\r) , (1.10) 

indicating that the fiavor branes are localized (all together) at the position r*= 0. Similarly 
the T^u[brane] contains delta functions with support on the position of the fiavor branes. In 
principle, one will need to solve second order, nonlinear, partial differential equations. 

Instead of directly dealing with the above equations, we want to present here an argument 
to understand which parameter controls the size of the fiavor effects on the action and, 
therefore, on the solution. We remark that the reasoning below is qualitative and in particular 
we will just write a background for flat Dp-branes as considered for instance in [21]. This will 
be enough for understanding the scaling with the parameters, at least in the cases studied 
in this review. In the following, we just focus on the behaviour with respect to Nf, N^, Qym 
and do not care about numerical prefactors. We will use notation similar to [21j. Consider 
the background associated to a stack of Dp color branes (in Einstein frame): 



7 2 I 

as = e 



a'U~P 



\ — 
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where Cp is a known numerical constant and U is an energy scale. On a background of this 
kind, we want to introduce Nj Dfc-fiavor branes and to know which is the relative importance 
of the associated terms in the action fll.6p and equations of motion (11.91) . With that aim, let 
us start by computing the coefficient in front of the term coming from the RR-form sourced 
by the color branes in (11. 6p . namely {2kIq)~^ y/gioe^'f' Fg_p. Using (II. 8p and (II. lip , we find 
that the lagrangian density associated to the color branes goes as: 

C^oior - {2k%)-'^o e^(a'-i)(^.iVe)(^-^)/2(v^t/)(^-3)(^-«)/2 . (1.12) 

Let us now look at the DBI term. We assume that the fiavor D/c-branes are extended along 
the Minkowski directions, the radial direction U and k — p — 1 directions within the sphere. 
We find: 

Cfiavor,DBi ~ NfTke~'''^/gk+i ~ — g^jf Ccoior , (1-13) 

where in order to get the last expression we have used (11.70 . (II. lip . (I1.12p and defined a 
dimensionless effective coupling as in [2T] : 

g'eff ~ 9ymNcU^-' ~ 9sN,{V^Ur-' . (1.14) 

Thus, parametrically, the action from the fiavor branes as compared with that from the color 

background is weighed by g^jj . We now want to take a low energy decoupling limit as 
in [21] (see also |22]), namely the dimensionless effective coupling g^ff and U are fixed as 
a' — )■ 0. Thus, the Veneziano scaling limit in this framework amounts to: 

Nc,Nf ^oo, geff fixed, Jf 3eff S^ed, (1.15) 

where the last relation comes from demanding that the fiavor effects are also fixed. Staying 
in the supergravity limit requires: 

1 « glff « iVj^ , (1.16) 

a constraint that limits the range of energy scales U for which the supergravity description 
is valid |2T]. Notice that if we further require that the fiavor terms do not parametrically 
dominate over the color ones, this can further restrict [/, depending on p and k. 

The probe limit, in which the fiavor action is negligible as cornpared to the gravity action, 
comes from making the last quantity in (ll.lSp vanishingly smalo As expected, that term is 
strictly zero in the 't Hooft limit. We now comment on the values of p, k that will appear 
in the following sections: 

•^In the literature, it is usually written that the probe limit is good when Nf ^ Nc- That is not strictly 
correct. For instance in the D3-D7 case, the probe approximation is valid when, parametrically, XNf <ti Nc- 
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For the D3-D7 case of section [21 the parameter weighing the flavor effects is {OYM-^cjjf ~ 
A^. For the cascading case of section the result is similar but one has to replace by 
the number of fractional branes M . Getting ahead of the discussion of upcoming sections 
we notice that, in these cases, it is not enough to take this parameter fixed, but it should 
also be small. This will be due to positive beta functions, as will be thoroughly discussed. 

From fll.lSp . we see that k = p (sections [3H and [6]) is particularly interesting^, since 
it is really jf- what has to be taken fixed. For this reason, only in these cases one can 
hope to describe - within gravity - phenomena as Seiberg-like dualities (see section l3.6.2p . 
Loosely speaking, the Klebanov-Strassler duality cascade [12] lies in this class of k = p 
theories, since it can be understood as the interplay between two sets of D5-branes wrapping 
vanishing two-cycles. 

We close this section by mentioning other brane intersections that will not be discussed 
further in later sections. In a D2-D6 system, the effective coupling f ll.l4p decreases at large 
U and therefore the flavor backreaction on the glue fades away in the UV - see ( ll.lSp - as 
expected in a superrenormalizable theory. This was observed in |23] when studying a solution 
with localized D6-branes. In a D4-D8 intersection, the opposite happens. The probe brane 
approach can be valid in an intermediate regime 1 ^ g^jj: <^ ^ but at a given U the 
fundamentals eventually take over and dominate. Notice that the value of U for which the 
D4-D8 theory loses its validity is parametrically smaller than that for which the unflavored 
D4-brane theory becomes pathological, which is set by f ll.l6p . 

1.4 The method 

Looking back at eq. (11. 6p . one can appreciate that in general finding the solution describing 
the backreaction between the type II closed strings and the open strings described by the 
Born-Infeld action is quite a challenging problem. Indeed, the fact that the flavor branes 
(BIWZ) are localized in the ten-dimensional space implies that we will have to solve second 
order, nonlinear, coupled and partial differential equations with localized sources. Basically 
what makes things so difficult are the presence of delta function sources and the fact that 
the differential equations describing the dynamics are "partial" (in principle depending on 
all the variables describing the space transverse to the flavor branes). In order to get some 
intuition of the answer, we may consider the case in which we will "erase" the dependence on 
these transverse coordinates (this is like considering the "s-wave" of the putative multipole 
decomposition of the full solution in this transverse space) and delocalize the sources. To 
achieve this, we will propose to smear the flavor branes over their perpendicular space. 

On the field theory side, this will amount to considering systems where the addition of 
the degrees of freedom transforming in the fundamental does not break any of the global 

Even if in all these sections we will deal with wrapped branes and therefore the backgrounds are not 
that similar to the argument above still yields the correct result. 
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symmetries of the unflavored QFT. Also, it may happen that the original U{Nf) is explicitly 
broken to U{l)^f as we are separating the flavor branes - see the discussion in section 7 
of [24j and in section 2 of [25]. An intuitive understanding of the smearing procedure will 
be discussed in section 11.51 while a more formal approach will be treated in section [7l For 
technical reasons, this procedure is cleaner in examples preserving some amount of SUSY, 
since the force between flavor branes is cancelled and the smearing is at no cost of energy. 
In the examples described in the following sections, we will proceed like this: 

• Consider an unflavored string background and find the embedding of flavor branes 
that will preserve some SUSY, or (in non-SUSY examples) that will be stable and 
solve the equations of motion for the brane. In the SUSY cases, this can be achieved 
by considering kappa-symmetric embeddings, that we review generically below. 

• Consider now Nf flavor branes in that particular embedding and smear them, getting 
an action in ten dimensions, as will be explained with all generality in section [71 

• Solve the equations derived from fll.6|) . that will now contain smeared branes and will 
be ordinary differential equations. In SUSY cases, there will be a set of BPS equations 
to be solved. In non-SUSY examples one might manage to get a fake superpotential 
and fake-BPS equations [26] . 

Moreover, one has to check that the flavor embeddings considered are still a solution in the 
backreacted geometry. 

Let us review briefly the main technical points collected above. 

1.4.1 BPS equations, Kappa symmetry (SUSY probes) and smearing 

Let us consider the case of a supersymmetric background, namely a solution of type II 
sugra for which the supersymmetry variations of the gravitino and the dilatino vanish 6ipn = 
6X = 0. We will not give here details on the form of these expressions, which can be found 
elsewhere. For instance, the string frame SUSY transformations of both type IIA and type 
IIB are written down in appendix A of [27] . 

Given a background that preserves some amount of SUSY, the idea is to find the hyper- 
surfaces in which to extend the flavor branes (in other words, finding the embeddings for 
flavor branes) so that these will preserve all (or a fraction) of the SUSY of the background. 

One then writes an eigenvalue problem, imposing that the preserved spinors of the back- 
ground are eigenspinors of the kappa-symmetry matrix 



= e . 



(1.17) 



See the papers 




for the definition of F 
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Once we have the kappa-symmetric embeddings as described above, we now proceed to 
write an action describing the dynamics of the closed and open strings, as in eq. (11. 6p . We 
then reahze that the problem will lead (unless we are adding D9-branes) to a system of partial 
differential equations. As discussed above, we then proceed to smear these flavor branes. 
For this we propose an ansatz for the metric, where the embedding of the flavor branes is 
clear and distribute them along the directions of their transverse space. This distribution 
of the flavor branes can be done in a uniform way. In some sense, we are 'deconstructing' 
the transverse space to the flavor branes by adding at each point one of the many Nf flavor 
branes. 

The key point is that once the BPS equations and kappa-symmetry conditions are simul- 
taneously satisfied, the problem is solved. In fact, it is a general result [29] that the SUSY 
equations Stpf^ = SX = 0, together with the Bianchi identities - and equations of motion - for 
the different forms modified by calibrated (namely, kappa-symmetric) sources imply the full 
set of equations of motion. 

In the following, we will discuss first an intuitive way of understanding this smearing. 
Then we will apply this to different examples in sections [2ll6l Finally, in section [71 we will 
present a formal way of implementing the backreaction from smeared sources. 

1.5 A heuristic viewpoint 

In the following sections, we will introduce the necessary mathematical machinery to con- 
sistently compute solutions of string theory in which smeared backreacting flavor branes are 
present. Before that, it is worth to make a digression in order to explain the general set-up 
in simple, heuristic terms. 

Let us make an analogy with electrostatics. Suppose we want to compute the electric 
field generated by a point-like charge and a couple of lines of charge, as depicted on the 
left of Figure [21 In order to depict the situation, we show dimension 1 lines of charge in 
a total space of dimension 2, but clearly the situation can be generalized by changing such 
dimensions. Since in the left plot there is no particular symmetry in the configuration, the 
resulting electric field will have a not so simple expression. But let us imagine that we 
consider a huge number of lines of charge as in the plot of the right, and homogeneously 
distribute them in the angle they form with the horizontal axis. In the limit of many lines, 
radial symmetry is recovered, the charge density is "smeared" and will be just given by a 
single (monotonically decreasing) function p{r). The electric field, accordingly, will also be 
radially symmetric. Notice that this process of describing a large number of discrete objects 
by a continuous distribution is ubiquitous in physics: for instance a "homogeneous" gas 
is a collection of atoms or the "homogeneous" Universe considered in cosmological models 
contains a collection of galaxy clusters. Also, solutions with different kinds of smeared 
sources have been considered many times in string theory contexts not necessarily related to 
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gauge-gravity duality, see for instance [30], [3T] . 




Figure 2: On the left, a point-like charge (in red) and two lines of charge at different angles. 
On the right, a configuration with many lines of charge. In the asymptotic limit of an 
infinite number of lines, they just correspond to a radial charge density. This picture depicts 
an analogous situation to the case of smeared flavored branes, when the fundamental fields 
are massless. 

When comparing to the string theory set-up, the point-like charge in the center corresponds 
to the color branes and the lines to the flavor branes (which have to extend to infinity) . The 
limiting, radially symmetric configuration, corresponds to the kind of smeared situations that 
we will discuss in this reviewj^l All functions of the ansatz can then be considered to depend 
on a single radial variable. For flavor branes, the different "angles" correspond to adding 
fundamental matter which couples differently to the rest of the fields. In some of the cases 
discussed in the following, we will see how this is reflected in the field theory superpotential 
(section EXH). 

We can still get further intuition from this simple analogy. In the right plot of figure |21 we 
see that all lines intersect at the center. From the string point of view, that means that the 
flavor branes are stretched down to the bottom of the geometry and the quarks are massless. 
In this situation, the charge density p(r) is highly peaked at r = 0. Essentially, that is the 
reason why for the solutions with massless quarks described in the following sections there 
is a curvature singularity at the origin, where all flavor branes meet. 

Then, a simple way of getting rid of such a singularity is to displace the lines of charge from 
the origin, while still keeping the radial symmetry. This is depicted in figure [31 If we dub the 
distance from any of the lines to the center as r^, the density of charge p{r) will vanish for 
< r < Tq, while it will asymptote to the "massless" = one as r ^ r^. From the brane 
construction, this displacement typically corresponds to giving a mass to the fundamentals 

^More precisely, it corresponds to the situations analysed in sections [2][5l For the cohomogeneity 2 cases 
analysed in section IH the different functions depend on two radial variables. A heuristic picture for such 
situations is presented in section [B] 
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(or, in particular cases, it could correspond to a non-trivial vacuum expectation value). The 
solution of section 12.3.21 which indeed is regular in the IR, is a neat example of this notion. 
Another possibility is to add temperature and to hide the singularity behind a horizon, see 
section 12.51 



Figure 3: These pictures depict analogous situations to the case of flavored branes, when the 
fundamental fields are massive. On the right, again, we add a large number of lines, such 
that in the limit radial symmetry is recovered. 

Going back to electrostatics for the right plot of figure [3l we know from Gauss' law that 
the charge density outside does not affect the central region. The corresponding statement 
in the field theory is that the massive fields decouple from the IR physics below the scale 
given by their mass. We find it interesting that, through this heuristic reasoning. Gauss' law 
is connected to the decoupling of heavy particles (or holomorphic decoupling in the SUSY 
cases) . 

Even if the example of electrostatics is useful to qualitatively picture what we will do in 
the following, the analogy is by no means perfect. We note two differences: first, we will 
be working with gravity, which is non-linear and, thus, one cannot find the final solution 
by superposing the fields generated by different sources (which in the case of electrostatics 
would make it rather trivial to find the electric field for the configurations depicted on the 
left of figures [2] and [3]) . Second, our "lines of charge" (the flavor branes) are dynamical. This 
means that is it not enough to compute the background fields generated from the sources 
but one also has to check that the sources are stably embedded in the geometry. 

We end this section by summarizing the pros and cons of looking for duals of unquenched 
theories for which the string solutions include smeared flavor branes, many of which can be 
inferred from the heuristic discussion above. On the positive side: 

• The smearing simplifies the situation allowing us to write ansatze depending on a 
single radial coordinate, and therefore the problem is eventually reduced to a set of 
ODEs. (For the cases of section [6l they depend on two radial coordinates and thus one 




16 



finds PDEs in terms of two independent variables, but again, without delta-function 
localized sources.) 

• Possible issues related to singularities and strong coupling are ameliorated in the same 
sense as they are washed out in electrostatics when considering a smooth charge density 
rather than a sum of delta-functions over a large number of electrons. 

• It allows a simple application of the powerful mathematical tools of calibrated geometry 
[32] . see section [71 

On the negative side: 

• Obviously, if we require the fiavor branes to be smeared, we are limiting ourselves to 
considering a very particular subset of all the possible fiavored theories. In particular, 
we require the superpotentials to effectively recover (some of) the global symmetries 
of the theory without fiavors. 

• Related to the previous point, one cannot realize, in general, theories with U{Nf) 
flavor groups. Since the flavor branes are required to sit at different points in the 
internal space, the typical string connecting different flavor branes is heavy and the 
flavor symmetry is typically broken to U{l)^f (one may also interpret the solutions as 
having flavor sjTumetry U{k)^i^^ for some k <^ Nf). From the point of view of the fleld 
theory, this amounts to having a Veneziano expansion with "one window graphs" , as 
pointed out in [25] . In principle, this fact can hinder the realization of some interesting 
physical features in the dual set-ups considered. 

1.6 Localized sources and other approaches 

As already remarked, this review focuses on solutions of string theory for which there are D- 
brane sources homogeneously smeared over a given family of possible embeddings, and that 
can be interpreted as duals of strongly coupled gauge theories in the Veneziano limit. As 
stressed above, this is a very particular subset of all the possible duals of theories with flavor. 
In a generic case, one should consider the sources to be localized at certain positions, such 
that the density of charge is given by a sum over Dirac delta functions. Such generic case is 
technically more challenging. However, remarkable works along these lines have appeared, 
pursuing solutions with the flavor branes localized at a single point of space (notice this is 
not the most general case either). We will not review them in any detail here, but the goal 
of this section is to provide a brief guide to the literature on the subject. 

The main ingredient of this approach consists of flnding solutions of supergravity which 
can be interpreted as intersections of branes of different dimension, with each stack of branes 
localized at a fixed position of space-time. In the context of gauge-gravity duality, the search 
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for such solutions was initiated in [33], [31] • These papers discussed D3-D7 intersections, 
which have been the most studied in the hterature (see below for different set-ups). A lot of 
progress was reported in subsequent years [35], [36], [37]. Among other aspects, these papers 
presented a better understanding of the D3-D7 solutions, the inclusion of fractional branes 
and clear matchings with field theory issues such as the running of couplings and anomalies. 
Further work on the D3-D7 localized system was performed in [38] (where the conifold was 
also addressed), [3S] (where D7 brane backreaction on bubbling geometries was considered) 
and [10], where the solution corresponding to D3-D7 in flat space was completed by providing 
an expression for the warp factor in closed form. It is also worth mentioning [JT] where a 
flavor D7-brane in a cascading theory was considered and its backreaction introduced as a 
perturbation. The finite temperature generalization of the background of [5T| was discussed 
in [12]. 

Let us now outline the literature on D2-D6 localized intersections, which can be interpreted 
as duals of 2+1 supersymmetric gauge theories coupled to fundamentals introduced by the 
D6-branes. The construction of the type IIA solutions (and their relation to M-theory) was 
carried out in [13], [B], [23]. In [15], meson excitations of this background were discussed 
and, in particular, the holographic dictionary relating meson-like operators to certain (closed 
string) supergravity modes was presented. On the other hand, the authors of [16] found a 
finite temperature version of the solution, which was used to discuss the thermodynamics 
of the system. Very recent progress in the D2-D6 systems, their M-theory uplifts and the 
detailed relation to Chern-Simons theories with flavor has been reported in [H]. 

Regarding D4-D8 intersections, localized solutions in that set-up were constructed in [48j 
in an early attempt to build a QCD dual. In the context of the Sakai-Sugimoto model [19] . 
backreaction from localized D8-D8 branes was analysed in |49j . 

It is also worth mentioning recent solutions in heterotic string theory which were argued 
to be related to flavored theories [50] . 

Interestingly, there are a few papers in which similar situations were considered in sub- 
critical string theory and therefore defined in dimensions lower than ten. In many of these 
cases, each flavor brane fills the whole space-time (therefore they are not localized, neither 
smeared). Some physics can then be extracted by using exact string theory methods but 
what these models have in common is that it is not possible to handle them within a well- 
controlled gravity description: gravity-like actions with just two derivatives suffer curvature 
corrections which cannot be neglected, nor consistently computed. However, there is the 
hope that the two- derivative actions can nevertheless provide additional non-trivial insights 
in the physics of the system. This idea was put forward by Klebanov and Maldacena in |51j . 
who considered a D3-D5 system in a six dimensional background (the cigar). Such a system 
is dual to 4D J\f = 1 SQCD as was shown using exact worldsheet methods in [52l |53]. For 
a recent discussion on the dual to the flavor singlet sector of A/" = 2 superconformal QCD 
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in a subcritical string framework, see [51]. The set-up of [5T] was generalized to different 
situations in [55], [TD], [56]. The finite temperature physics of a model in [10] was analysed 
in [57]. Bottom- up approaches (in which a high-dimensional gravity theory is proposed to 
describe some specific features of QCD) with space-time filling flavor branes have been dis- 
cussed in [HE], [59]. Recently, a bottom- up approach to the conformal window along these 
lines has appeared [60] . 

Finally, let us mention a recent contribution by Armoni [61] , in which a way of departing 
from the quenched approximation was proposed. The fermion determinant is expanded in 
terms of Wilson loops. It then turns out that a sum of correlators of an observable with 
the Wilson loops boils down to an expansion in which can in principle be computed. 
It would be nice to further develop possible implications of this observation in holographic 
set-ups. 

2 Flavor deformations of AdS^ x 

Our first concrete application of the procedure described will be the flavor deformation of 
74^5*5 X 5*^. This is the simplest possible case and, hopefully, it will neatly illustrate the 
comments of section 11.21 In fact, for most of the discussion, the formalism applies to any 
AdS^ X geometry, being a five-dimensional compact Sasaki-Einstein (SE) space, so 
we will refer to this more general case during this whole section. At some points, we will use 
the two notable examples X^ = or Ti^i to clarify particular issues. 

Let us start with a general comment. Since the AdS^ x theories without flavor are 
conformal, we expect that once we include extra matter, a positive beta function is generated. 
This is in fact the case and leads to the appearance of a Landau pole. Nevertheless, as in 
QED, the theory renders meaningful IR physics even if the UV is ill-defined, as long as the 
IR and UV are well separated scales. However, this separation does not allow to have Nf 
and Nc of the same order. As we will see, one can define a parameter e ~ XNf/Nc which 
weighs the internal flavor loops and that has to be kept small. The effect of the unquenched 
quarks can then be computed as an expansion in e. 



After introducing the framework in sections 12. H 12. 2[ we present the unquenched super- 
symmetric (A/" = 1 in 4d) solutions in 12.31 In section \2A\ we present an instance of the 
effects of the unquenched flavors on a physical quantity, namely on the mass of a particular 
meson tower. Then, in section 12. 5[ we break supersymmetry by turning on temperature and 
analyse the physics of the dual quark-gluon plasma. We end in section [2l6] by discussing the 
range of validity for the solutions and approximations used. 
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2.1 The geometries and field theories without flavors 



The models we discuss here are obtained by placing a stack of A^^^ D3-branes at the origin of 
the six-dimensional cone over X^. The corresponding type IIB background reads: 

F5 = dh^^ dx^ A ■ ■ • A dx^ + Hodge dual , 

h(r) = ^ Q = (21) 

4r4 ' - Vol(X5) ' ^^-'^ 

where we have taken the near horizon limit. The dilaton is constant and all the other fields 
of type IIB supergravity vanish. In general the metric of the SE space can be written as 
a Hopf fibration over a four- dimensional Kahler-Einstein (KE) manifold: 

ds\s = ds^E + {dr + AkeY , (2.2) 

where r is the fiber and A^e is the connection one-form whose exterior derivative gives the 
Kahler form Jke of the KE base: 

dAKE = 2Jke . (2.3) 

Let us first consider the particular case in which X^ is the five-sphere S^. In this case 
the KE base is the manifold CP^ (with the Fubini-Study metric) and the space transverse 
to the color branes, with metric dr"^ + c?s|5, is just M^. When X^ = the coefficient Qc 
appearing in (12. ip is just Qc = IGirgsa'"^ Nc- Moreover, as is well-known, the field theory 
dual to the AdS^ x 5*^ background is A/" = 4 SYM in 4d, which, in M = 1 language, can 
be written in terms of a vector multiplet and of three chiral superfields {i = 1,2,3) 
transforming in the adjoint representation of the gauge group and interacting by means of 
the cubic superpotential: 



WV=4 = Tr 



$l[$2,$3] • (2.4) 



If we represent the transverse of the AdS^ x solution in terms of three complex variables 
Zi {i = 1, 2, 3), one can regard the Zj's as the geometric realization of the adjoint superfields 

The second prominent example which we will analyze in detail is the one in which X^ is 
the T^'^ space with metric: 



ds 
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where the range of the angles is ip E [0,47r), ipi G [0,27r), 9i G [0,7r]. Since Vol(T^'^) = 
167r^/27, the coefficient Qc for this solution is Qc = 27TTgsa''^ Nc. In this case the space 
transverse to the color branes is the conifold, which is a 6d Calabi-Yau cone which can also 
be described as the locus of the solutions of the algebraic equation: 

Z1Z2 = z^Zi , (2.6) 

where the Zi are four complex coordinates. The relation between these variables and the 
coordinates used in (12. 5p is the following: 

^-,=r2e2^^ sm — sm — , ^2 = ^^^ea^^^^^^^^^ cos — cos — , 

^ 2 2 2 2 

^3 = ef (^+^1-^2) cos y sin ^ , z^ = rf e5(^-^^+^2) y cos ^ . (2.7) 

Notice also that the metric written in (12. 5p is of the form (12. 2 p where the KE base is just 
the S"^ X space parameterized by the angles {6i, ipi) and one should make the following 
identifications: 

r = V'/S , ^Ti'i = o (cos0i(i(/?i + cos62d(p2) , 

J^i,i = ^"^J''' = -i ( sineidOiAdipi + sin 62 d92 A dip 2) ■ (2.8) 
2 6 

The field theory dual to the AdS^ x T^'^ background is the JV = 1 superconformal quiver 
gauge theory with gauge group SU{Nc) x SU{Nc) and bifundamental matter fields Ai,A2 
and -81,-82 transforming respectively in the {Nc,Nc) and in the {Nc,Nc) representations of 
the gauge group ^-e. the so-called Klebanov-Witten (KW) model. The matter fields 
form two SU (2) doublets and interact through a quartic superpotential 

Wkw = h e'^e"' Tr [ABkA.Bi] . (2.9) 

The fields Ai and Bi can be related to the coordinates Zi by means of the following relations: 

zi = AiBi, Z2 = A2B2, Z3 = AiB2, Z4 = A2Bi, (2.10) 

which automatically solve the defining conifold equation (12.60 . 



2.2 Flavor branes and smeared charge distribution 

The flavor branes for the AdS^ x backgrounds just described are D7-branes extended 
along the four Minkowski directions as well as along a non-compact submanifold of the cone 
over X^. The type of flavor that the D7-branes add depends both on the space X^ and on 



21 



the submanifold they wrap in the transverse space. We first illustrate the situation with the 
two examples of = and Ti_i and at the end display the general expressions. 

The first instance is the case in which = S^. In this case a simple kappa symmetry 
analysis shows that, in order to preserve eight supersymmetries, the D7-branes must be ex- 
tended along a codimension two hyperplane in which, in terms of the complex coordinates 
can be written as: 

aiZ^ + a2Z^ + asZ^ = fx , (2.11) 

with the Oj and /x being complex constants satisfying Yli = 1- the field theory side 
these flavor branes introduce M = 2 fundamental hypermultiplets (Q^, Qr) ('^ = 1? ■ ■ ■ ^/) - 
nonetheless, a generic collection of branes within the family (12. lip retains just M = 1 susy. 
The corresponding superpotential for an embedding such as the one in (12. lip can be written 
as: 

W = WV=4 + Qr [ 5^ a,- + m] Q^ (2.12) 

j 

where the mass m is related to the constant in (12. lip . Notice that since the embeddings 
are holomorphic, it is not possible to smear them in a way in which the full 5*0(6) isometry is 
realised. After smearing over the embeddings (12. lip , one can recover, at most, SU{3) xU{l), 
as will be seen directly from the dual solution. 

In the case of the AdS^ x T^'^ background there are two classes of holomorphic embeddings 
which correspond to different types of flavors in the KW theory. In terms of the Zi coordinates 
of (12. 7p the representative embedding of the first class is given by the equation zi = /i. This 
is the so-called Ouyang embedding [3Tj, which has two branches in the massless limit /i = 0. 
In each of these branches the D7-brane adds fundamental matter to one of the two nodes of 
the KW quiver and antifundamental matter to the second. The corresponding superpotential 
contains cubic couplings between the quark fields and (ji {i = 1, 2) and the bifundamental 
fields Ai and Bi. For example, for the massless embedding zi = the superpotential (12. 9 p is 
modified as: 

W,,=o = Wkw + hqiAiq2 + h2q2Biqi , (2.13) 

where, here and in the following, traces over color indices and sums over the Nf flavor indices 
are implied. The second class of D7-brane embeddings is the one giving rise to non-chiral 
flavors, whose representative element is given by the equation zi — Z2 = yU. In this case 
every D7-brane adds fundamental and antifundamental flavor to one node of the KW quiver 
and the flavor mass terms do not break the classical symmetry of the massless theory. The 
corresponding superpotential contains only mass terms and quartic couplings, namely: 

W = Wkw + h - Aa^ajgi + h2 U^iAi - ^s^jgs + {q.qif + m {q,qi) . (2.14) 

In order to develop our program and construct backreacted gravity solutions for smeared 
distributions of flavor branes following ref. [63], we should be able to find a family of 



22 



equivalent embeddings for each type of configuration described above. In the case of the 
AdS^ X background eq. (12. lip provides such a family. Notice that, even if each individual 
embedding of the form (12. lip preserves = 2, which supersymmetries are preserved depends 
on the Oj's. Nevertheless, one can check that all the holomorphic embeddings of the type 
(12. lip are mutually supersymmetric and, due to the holomorphic nature of the linear equation 
(12. lip , they preserve the same common four supersymmetries (A/" = 1) for all values of the 
constants Oj. Thus, we can use these constants to parameterize the family of different planes 
that constitute our continuous distribution of flavor branes. 

In the case of the AdS^ x T^'^ background one can generalize the chiral embedding zi = fi 
by acting with the SU{2) x SU (2) symmetry of the conifold. The corresponding family of 
embeddings takes the form: 

4 

^aiZi = n, (2.15) 
1=1 

with the complex constants a, spanning a conifold (up to overall complex rescalings) 

q;iQ;2 — «3«4 = . (2-16) 

Notice that embeddings like zi — Z2 = are not in this family. Indeed, the non-chiral 
embeddings zi — Z2 = fJ^ can be generalized as: 

pzi - pz2 + qz3 + qz4 = fi, (2.17) 

where p,q span a unit 3-sphere, i.e. they satisfy |pp + |gp = 1. 

In spite of the differences among the cases presented above, the charge distribution gener- 
ated by these families of embeddings can be written in a common form. The reason for this 
universality is the underlying Sasaki-Einstein structure. In order to illustrate this fact, let 
us consider the chiral embeddings (I2.15P in the case in which the mass parameter /i is zero. 
Without loss of generality we can rescale the coefficients and fix ai = 1. Then, (12.160 
fixes a2 = 0:304 and, after using (12. 7p . the massless embedding equation 

zi + asa4Z2 + asZ3 + a4Z4 = , (2.18) 

nicely factorizes as: 

(6 6 \ f 6 6 \ 

sin ^ + oge^'^i ^ ) ( ^ + "^e*'^' ^ ) " ° ' ^'^'^^^ 

Notice that the vanishing of each of the factors in (I2.19P determines a branch in which the 
branes sit at a fixed point of one of the two two-spheres parameterized by the angles (^j, (/9j). 
The constants 03 and 04 determine the particular point at which each brane is sitting in 
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each S"^. Indeed, if and ^2 systems of worldvolume coordinates for the D7-branes, 
these two branches can be written as: 

$,1 = {x^,x^,x'^,x^,r,92,'f2,^} , Oi = 91 = const. , yji = v?^ = const., 

= ^""^5 ^'^5 V^i5 ^} 5 02 = 02 = const. , (^2 = ^2 = const. (2.20) 
In figure H] we have represented the two branches for the embedding fl2.18l) . From the 




Figure 4: We see on the left side the two stacks of Nj fiavor-branes localized on each of their 
respective S^'s (they wrap the other S*^). The flavor group is clearly U{Nf) x U{Nf). After 
the smearing on the right side of the flgure, this global symmetry is broken to U{1)^^~^ x 

f/(l)^/-l X f/(l)B X Uil)A. 

fleld theory side, which particular embedding we choose determines the coupling between 
the associated quarks and the bifundamentals. Roughly speaking, the contribution to the 
superpotential of an embedding determined by some 03, 04 is /iigi(Ai + 04^2)52 + ^2Q'2 (-Si + 
a3B2)qi. Thus, when we smear and sum over all the possible and a^, both S'[/(2)'s (the 
one rotating the Aj's and the one rotating the -Bj's) is effectively recovered. Figure H] is the 
geometric interpretation of this effect. 

It is straightforward to compute the charge density produced by this localized D7-brane 
configuration. Indeed, taking into account the contribution of the two branches, one gets: 

Qloc ^ ^(2) _Q*^^^_ /\ ^ ^(2) /\ _ (2.21) 

To produce an homogeneous configuration of Nf D7-branes we should distribute in every 
branch the branes homogeneously along their transverse two- sphere. In the continuum limit 
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Nf — )■ oo this procedure amounts to performing an integration over each S"^ with the corre- 
sponding volume element, namely: 



n 



47r 



6*1, (fii- (fil) del dfi dOi A d(fi + 



+ 



^ sin e; 6^^^ [02 ~e;,<^2- <^;) de; d<^; 



d92 A d(p2 



(2.22) 



The integrations over 6* and (p* in f l2.22p can be immediately performed, yielding the fol- 
lowing expression for the smeared charge distribution of D7-branes: 



Nf ( 

= — - sin 6i dOi A du:>i + sin 62 d92 A du:>2 
An \ 



(2.23) 



Notice that in f l2.22p we have included the normalization factor Nj/An in such a way that the 
resulting distribution densities sin^* Nf/An are normalized to Nf when they are integrated 
over 5*2. Notice that, as already pointed out above, the flavor symmetry of the smeared 
configuration is f/(l)^-'' rather than U{Nf), since the branes are not placed on top of each 
other. Interestingly, a similar calculation for the embeddings fl2.17p in the massless case 
fi = gives rise to the same charge density for the smeared configuration as in ( ]2.23p [64j. 
This is because the form of Q in f l2.23p is determined by the SU{2) x SU{2) x Z2 global 
symmetry which we want to recover after smearing. 

Actually, one can rewrite fl2.23p in a form which can be easily generalized to any continuous 
family of equivalent D7-brane massless embeddings in an arbitrary Sasaki-Einstein manifold. 
Indeed, by using ( 12. Sp one can rewrite the right-hand side of ( ]2.23p in terms of the Kahler 
form of T^'^ as: 

?>Nf 



n 



In 



J' 



For an arbitrary Sasaki-Einstein space X^, the expression fl2.24p generalizes to: 

fi's ^ = ~2 Qf Jke , 
where Q/ is the following constant coefficient: 

Vo\{X^)gsNf 



Qf 



4Vol(X5) 



(2.24) 



(2.25) 



(2.26) 



In fl2.26p is the compact submanifold of wrapped by the D7-brane in a massless 
embedding . Notice that in this case the D7-brane worldvolume along the space transverse 
to the color branes is always of the form I x X^, where / is a non-compact interval along the 
holographic radial direction. It is worth noticing that the factor Vol (X^) /Vol (X'^) appearing 
on the right-hand side of fl2.26p is just the volume transverse to any individual flavor brane, 
over which we are distributing the D7-branes. For the massless chiral embeddings in the 
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conifold one can readily check, after taking into account the contribution of both branches 
in (ESO]), that Vo\{X^) = 1Qit^/9. Since Vo\{T^'^) = 167r727, one can easily prove that 
fl2.25p reduces to fl2.24p . In the case = the three-manifold X^ is just a unit 5*^ and 
Vol(X^) = 27r^. Therefore, we obtain the following values of Qf for X^ = S^,T^'^: 



f^^ V5 _ c-5 



2-K 



for = 5% 



Qf = { (2.27) 
M/ for = Ti'i . 

The charge density VL determines the ansatz of Fi in the backreacted geometry. Indeed, the 
WZ part of the D7-brane action contains a term in which the RR eight-form potential Cg 
is coupled to the D7-brane worldvolume. The continuous limit for this term amounts to 
performing the following substitution: 



Swz = TrJ2 Cs^Tr Q A Cs , (2.28) 



-J Ms, J Mio 

which leads to the following violation of the Bianchi identity for Fi. 

dFi = -2kIio) Ti^ = -9s^ ■ (2.29) 

Taking into account the general expression of ^2 for a massless embedding written in fl2.25p . 
as well as the relation (12. 3p between the one- form Ake and the Sasaki- Einstein Kahler form 
Jre, one is led to adopt [63] the following ansatz for F^. 

Fi = Qf{dr + AKE) . (2.30) 

A simple modification of ( 12.30p for Fi allows us to extend the ansatz to the case in which 
the quarks are massive [63j. This modification corresponds to introducing a function p{p) of 
the holographic coordinate p and performing the substitution Qf ^ Qfp{p) in fl2.30p . This 
function p{p) encodes the effects of the non-trivial profile of the D7-branes. Indeed, when the 
quarks are massive the brane does not extend along the full range of the radial coordinate 
p and, accordingly, p{p) must vanish for p < pg, where p = Pg is the radial location of the 
tip of the D7-brane. Moreover, the function p{p) should approach the value p = 1 when 
p ^ Pg since in this region the quarks are effectively massless. The form of the function p{p) 
is not universal and depends on the particular embedding of the D7-brane. For the three 
embeddings in the cases X^ = and T^'^ discussed above, the expressions for p{p) are given 
below. At this point let us simply notice that the charge density Q is modified with respect 
to the massless 

QsQ = -2 p{p) QfJKE - Qf Pip) dp A [dr + Ake) , (2.31) 
where the dot denotes derivative with respect to the radial variable p. 
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2.3 Backreacted ansatz and solution 



Let us now write an ansatz for the backreacted D3-D7 background for a generic Sasaki- 
Einstein space [63]. It is clear from the discussion of the previous subsection that, after 
performing the smearing, the resulting RR one-form Fi introduces a distinction between the 
directions of the U{1) fiber and of the KE base of X^. Therefore, it seems clear that the 
effect of the smeared flavor branes on the metric should be an internal deformation of the X^ 
in the form of a relative squashing between the KE space and the Hopf fibej^. Accordingly, 
let us adopt the following ansatz for the metric in Einstein frame: 



h{p) 1 ' dxl^ + h{p) ' e^f^P^dp" + e2s(p) ds],E + e^f^P^ {dr + A 



KE) 



(2.32) 



where g{p) and /(p) are the functions that implement the squashing mentioned above and 
the function multiplying dp^ amounts to choosing a particular radial variable p which is 
convenient for our purposes. Moreover, the dilaton will depend on p and the RR forms 
and Fi have the form: 



= 0(P) 



F5 = Qc(l + *)e{X' 



Fi = QfP{p) {dr + AKE) 



(2.33) 



where s{X^) is the volume element of X^ and Qc and Qf are written in (12. ip and (I2.26P 
respectively. The function p(p), whose form depends on the D7-brane embedding, takes into 
account the effects of massive quarks, as explained above. 

Given the ansatz (I2.32I) - (I2.33I) one can easily study the supersymmetric variations of the 
dilatino and gravitino in type IIB supergravity and find the corresponding first-order BPS 
equations, which ensure the preservation of four super symmetries. The resulting equations 
are [631: 



d,f = 3- 2e'f-^^ - ^p{p)e^ 



QfP{py 



dph 



-Qce 



-ig 



(2.34) 



Remarkably, the system (I2.34p can be integrated analytically for any function p{p). In order 
to present this solution, let us define the function 77 (p) as follows: 



(2.35) 



where Pq is the value of the radial coordinate at the tip of the flavor brane (p(p < pq) = 0). 



^Just in the case when is the sphere S^, this squashing breaks part of the isometry 5*0(6) — > SU{3) 
U{1), where SU{3) is the isometry of the Kahler-Einstein base CP^. 
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Then, we can write down quite simple expressions for /, g, 0, namely: 



e 



^ - cae^e-e (l + e-^''(cie^ + r/)) ^ , (2.36) 



where we have introduced a reference scale p^: and we have defined 0* = (f){p = p*). Notice 
that the warp factor h can be obtained as the integral of e~^^, as follows from the last 
equation in the BPS system fl2.34p . In fl2.36p Ci and Ca are integration constants that we 
now fix. First, if we demand IR regularity of the solution, we need g = f when p < pg. Since 
r] vanishes at p = p^, we need Ci = 0. Moreover, the constant Ca is just some overall scale 
and has no physical meaning. It is natural to fix it to a'^e^ in order to give appropriate 
dimensions and to recover the usual expression for the metric when Qj = 0. Therefore, we 
find: 



& = ya e^e 



^ (1 + e-^Pj]) ' . (2.37) 



Notice that, when = 0, we recover the unflavored AdS^ x background. Indeed, in this 
case = 0* and i] = and, after performing the change of the radial variable r = a/o'c'', 
we get that = = r and the background fl2.32p - fl2.33p coincides with the one written in 

Let us now introduce the following parameter 

= Qfe'^* , (2.38) 

which, as we will see in a while, controls the effects of quark loops in the backreacted 
supergravity solution. Indeed, the gauge/gravity dictionary for the type of theories we 
are studying relates the exponential of the dilaton to the Yang-Mills coupling constant. 
For example, for the (flavored) Af = 4 SU{Nc) theory, dual to the deformed AdS^ x 
background, the gauge coupling is gyM = 4:7i ggC'^ and, thus, the 't Hooft coupling at the 
scale p=K is given by: 

A, = 4:71 gsN^e'f'' . (2.39) 

For the quiver theories that correspond to different geometries, the gauge groups are of 
the form SU {Nc)"". Let us generalize a relation from the orbifold constructions Angy^n = 
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{9sG'^)~^ [621 165], and consider all the gauge couplings gYM,i to be equal. Then Air ggNcC'^, 
strictly speaking, gives the 't Hooft coupling at each node of the quiver, divided by n. 
However, with an abuse of language we will simply refer to it as the 't Hooft coupling. 
Therefore, by using ( ]2.39p and the definition of Qj in ( 12.26^ in ( 12.38p . we get: 

VoUX'^) Nf , , 

^ ' X.-f. (2.40) 



c 



16n Vol{X^) N, 
In particular, when = this relation becomes: 

e*(x^=5^) = ^A*^- (2-41) 

Notice that the fact that (f) is not constant in the backreacted solution is simply a reflection, 
in the gauge theory dual, of the running of the Yang-Mills coupling constant when matter 
is added to a conformal theory. 

In terms of the dilaton and the function rj of f l2.35p take the form: 

- ^T77V37 ' V = e*e<^-^* T e'^ piO ■ (2-42) 



1 - £ P(0 d^ 



p, ^^^^ ^ ^Pg 



One of the prominent features of our solution is the fact that, for Nf ^ 0, the dilaton blows 
up at some UV scale p = plp, determined by the condition: 



Plp 



piO dC = e;' . (2.43) 

Clearly, in order to have a well-defined solution, we should restrict the value of the radial 
coordinate p to the range p < pip. In view of the relation between the Yang-Mills coupling 
gYM and the dilaton {gyM ~ the divergence of implies that gvAi blows up at some UV 
scale, i.e. that the gauge theory develops a Landau pole. This UV pathology of our solution 
was expected on physical grounds since the flavored gauge theory has positive beta function. 
Indeed, we will check below in some particular case that our solution reproduces the running 
of the coupling constant of the dual field theory. 

2.3.1 The supersymmetric solution with massless quarks 

We now consider the particular case of massless quarks, which corresponds to taking the 
charge distribution given by f l2.25p or simply p{p) = 1. In this case fl2.42p simply gives: 

e^-f* = -— ^ . , e-'^v = 'ie^-'^* , (2.44) 

1 + (p* - p) 6 
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and the solution written in f l2.37p reduces to: 



= y/^ eP [L + e^[p^ - p))'^ \ L + e^[- + p^ - p 



+ + (2,45) 



dh Qc 
dp 

Notice that the location of the Landau pole in this case is just plp = p* + e^^ and that the 
range of p for which the solution (I2.45P makes sense is p G (— oo,plp). Moreover, by using 
the definition of 6=,, in fl2.38p one can immediately show that the dilaton can be written as: 

e'^(^) = r . (2.46) 

Qf{pLP-p) 

Let us now verify that the dependence on p of in f l2.46p matches the expectations from 
field theory. For concreteness we will consider the case of A/" = 4 SYM with matter. Similar 
checks can be done in other cases (see [63] for the case of the Klebanov-Witten theory). By 
using the relation between the Yang-Mills coupling and the dilaton discussed above, as well 
as the value of Qj for = written in fl2.27p . one gets: 

^ = Nf ipLP - p) . (2.47) 
9ym 

In order to read the running of the coupling constant from fl2.47p we must convert the 
dependence on the coordinate p in fl2.47p into a dependence on the energy scale of the 
corresponding dual field theory. At an energy scale p much lower that the Landau pole scale 
Ajjv (^-e. for p << p^p) the scaling dimensions of the adjoints and fundamentals take their 
canonical values and the natural radius/energy relation is: 

PLP- P = log ^ . (2.48) 
A* 

Plugging this relation in f l2.47p we get: 

^ = iV,log^. (2.49) 

9ym 

Therefore, we get a logarithmic scaling of the coupling of the type Sn^/gYM = blogE, 
with b = —Nf, which matches the one- loop field theory result^ in which one has that b = 

^In principle, one could object that, being strongly coupled, the matter fields could get large anomalous 
dimensions making this result suspicious. However, since we are performing a small perturbative (in e*) 
deformation of the unflavored backgrounds, the anomalous dimensions for the fundamental multiplets cannot 
differ much from their quenched values. For the = case, those anomalous dimensions vanish. We 
thank F. Bigazzi for stressing this point to us. 
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- 3N^ - Nf. 

In order to have a clearer understanding of the deformation of the AdS^ x metric 
introduced by the flavor, it is very convenient to change to a new radial variable r, which is 
defined by requiring that the warp factor takes the same form as in the unflavored case (see 
(O)): 

Qc 
4 

By integrating the last equation in fl2.45p we can get h{p) and thus r(p). We will perform 
this integration order by order in a series expansion in powers of e*. The additive integration 
constant will be fixed by requiring that r(p*) = r^, = y/a' e^* . One gets: 



(2.50) 



a '2e' 



5e 



10368 



-(19-24(p,-p) + 144(p, -p)^ 



6e4p-4p*(3 + 4(p,-p)) 
(2.51) 



It is now straightforward to obtain the functions f{r), g{r) and the dilaton 0(r) as expansions 



in powers of e^,. Up to second order we have: 

/I 



of 



r 



1 4 

— (1 H ) H — 

2V 3rr 1152 



17 



94 5 
+ 9^:8 



48 log( 



24 



Ir^ 



+ 



1152 



106 r'^ 5r^ , , r 



r e 
+ log - + 

r* 72 



1 - ^ + 12 log - + 36 log" - + 0(e 



(2.52) 



Eq. f l2.52p neatly displays the effects of quark loops in the deformation of the geometry and 
in the running of the dilaton (the latter is related to the running of the gauge coupling, as 
argued above). It is important to point out that the deformed geometry has a curvature 
singularity at the origin r = (or p = — oo) (this singularity is similar to the one that 
appears at r = in a 2-dimensional manifold with metric ds"^ = dr"^ + r^(l + r)dip'^). In 
the same IR limit, e"^ runs to 0. As argued in section [T31 the appearance of this singularity 
can be intuitively understood as due to the fact that, in this massless case, all branes of our 
smeared distribution pass through the origin and the charge density is highly peaked at that 
point. From the field theory side, one can think of the singularity as appearing because the 
theory becomes IR free, as first pointed out in [M]. Consistently with these interpretations 



and with the heuristic picture of section II. 5[ the IR singularity can be easily cured by giving 
a mass to the quarks (it is a "good" singularity according to the criteria of ^6], [67]). We 
will explicitly verify this fact in the next subsection. 
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2.3.2 The supersymmetric solution with massive quarks 



Let us now find the backreacted supergravity solution for massive quarks. As mentioned 
above, the function p{p) entering the ansatz for Fi in this case is not universal and depends 
on the particular Sasaki-Einstein space and on the family of D7-brane embeddings cho- 
sen. For concreteness we first concentrate in discussing the case in which = S^. The 
calculation of the function p{p) in this case was performed in appendix C of [68]. If = e''«, 
one has: 



I _ g2(pg-p) 



0(p-p,) 



(2.53) 



When p > Pq the function p{p) is non-vanishing and one has to perform the integrals ap- 
pearing in f l2.42p . These integrals can be straightforwardly done in analytic form and yield 
the result: 



1 



+ p* - p 



3 
— ( 
2 



p6p<,-6p _ ;:^g2p,-2p _| g 



1 



ipq-ip _ ^^Apg-ip, _|_ g2p,-2p. 



6 ' ' 6 



(l + e*(| + - p - |e6P9-6p _ |e2p9-2p + |e4p,-4p _ Ig4p5-4p* _^ g2p,-2p,)^ a 



1 ) 

3 



log 



1 + e* p=^ - p 



,2p,-2p 



+ 7^ 
4 



4p9-4p _|_ g2pg-2p, _ j!^^4pq-4p. 



(2 



As a check, notice that setting pg — > — oo one recovers the massless solution of fl2.44p and 
f l2.45p . We still have to write the solution for p < pg. In this case p{p) vanishes and the 
dilaton is constant and, by continuity, it has the value that can be read from fl2.54p inserting 
P = Pq- 



4>IR = 0g = 0* - log(l + (p* - Pq 

The functions / and g are equal and given by: 



_ _L p2p,-2p, 

4^ 



4 



4pg-4p, 



)) 



a 2 e'^e « 



(P < Pq) ■ 



(2.55) 
(2.56) 

It follows straightforwardly from these results that the IR singularity at p = — oo of the 
massless case disappears when p ^ since the background reduces to AdS^ x for p < pq. 
Moreover, one can verify that the metric is also regular at p = pg. Thus, as stressed in 
section 11.51 the smearing of massive fiavors allows one to smooth out IR singularities. 

Similar calculations can be done for the conifold theories. In this case we redefine the 
parameter p of the embedding equations (I2.15P and (I2.17P as |p| = e~ . The charge distri- 
bution for the family fl2.15p of chiral embeddings was obtained in ref. [25], with the result: 



p{p) = 1 - e'^P^-P^ (1 + 3p - 3p,) e(p - Pq) . 



(2.57) 
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Similarly, for the non-cliiral embeddings f l2.17p the function p{p) is given by [69] : 



P(P) 



I _ g3(pg-p) 



0(p-p,) 



The corresponding supergravity solutions have been written down in refs. 
regular in the IR, much in the same way as in the = case detailed above. 



(2.58) 
They are 



2.4 Screening effects on the meson spectrum 

The holographic theories with flavors present mesonic excitations, meaning that there exists 
a spectrum of colorless physical states created by operators which are bilinears in the fun- 
damental fields. They are associated to normalizable excitations of the flavor branes as was 
neatly explained in the seminal paper [6j. For a review of this broad subject, see |7j. Notice 
that the notion of "meson" we use here generalizes that used in QCD. For instance, the 
"mesons" of [6] are excitations of a non-confining theory and in this case the dimensionful 
quantity that sets the meson masses is just the quark mass (divided by a power of the 't 
Hooft coupling), not a dynamically generated scale. 

In the present section, we review how the presence of unquenched flavors can affect the 
discrete mesonic spectrum. Again, we will restrict ourselves to the smeared set-up and follow 
[70] . For discussions about screening effects on the spectrum in cases with localized rather 
than smeared backreacting flavor branes, we refer the reader to [lOl [7]. The effect of the 
smeared flavors on the hydrodynamical transport coefficients (in a finite temperature setting) 
was studied in [681171] . It is also worth mentioning that, within the model we will introduce 
in section [3l the screening effects on the glueball spectrum have been recently analyzed in 

For the sake of briefness, we will just focus on an example and discuss a particular mesonic 
excitation in the backreacted Klebanov-Witten model. The analysis and conclusions for 
different modes and/or different models should be similar, see [70] for some other examples. 
In particular, we will consider oscillations of a D7-brane which introduces massive non-chiral 
flavor [73], and just look at the oscillation of the gauge field that gives rise to a vector mode 
in the dual field theory. Thus, we discuss the physics of a meson whose "constituent quarks" 
are massive in the presence of many dynamical massless flavors. 

We write the gauge field along the Minkowski directions as = a^^p) C,^e^^^ , where 
is a constant transverse vector. The equation describing this oscillation was found in |70j . 
building on the formalism introduced in [73]. It reads: 

= dp (e23-3^(e3^ - e^P^)d,a,) + M^he^^+^-f (^1 + gSPQ-sp ^^2^-2/ _ ^ (2.59) 

where = —k'^, the constant pq is the minimal value of p reached by the D7-brane (related 
to the quark mass) and f,g,h are given in f l2.45p . 
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Notice that for the meson excitation, we just use a D-brane probe, namely, we consider the 
oscillation of a single brane in a fixed background. At first sight, this could look contradictory, 
since our aim is always to take into account the effect of the flavor branes on the geometry. 
Then, one may think about considering coupled fluctuations of brane and background fields. 
Nevertheless, this is not necessary: there are Nf 3> 1 flavor branes which are affecting the 
background but when we consider a meson, only one (or two) out of this Nf is fluctuating. 
Therefore, the effect of this oscillation on the background is suppressed by A^^^^ with respect 
to the contribution of the whole set of branes and therefore is consistently negligible. On 
the other hand, the existence of the rest of flavors - and the associated quantum effects on 
the spectrum - are taken into account through the deformation they have produced in the 
background geometry. 

Following the standard procedure [HI [7] , a discrete tower of values for should be found 
when selecting solutions of (12.591) which are regular and normalizable. Since the background 
has a Landau pole, some prescription is needed for dealing with the UV limit (large p). 
Technically, we will just require that the fluctuation vanishes at p^,. Physically, one 
can check that this is a consistent procedure if pq « p*: we are interested in some IR 
physics which should be independent of the UV completion of the theory at p > p*, up to 
corrections suppressed by powers of the UV scale. Namely, we neglect contributions of order 
qPq-p* ^ Am. and check that the spectrum can be written in terms of IR quantities. The 



value p^ disappears from the final result, apart from the quoted negligible corrections. See 
[70] for further discussions on the issue. In section 12.51 we will see similar examples of how 
to deal with the Landau pole. In that case, the IR scale, which has to be much smaller that 
the arbitrary UV scale at p^,, is set by the temperature rather than by the quark mass. 

In order to estimate the spectrum from (12.591) . we can use a WKB approximation. In 
|70] . using a formalism developed in [71], an expression for the mass tower in terms of the 
principal quantum number n was found. Adapting notation to the one we are using here: 



Let us evaluate this integral at first order in e*, by inserting (I2.45p . We still have to fix 
the additive constant for h, which we can do by requiring /i(p*) = (in [70] h{pLp) = 
was used. It is crucial that both prescriptions give the same result, up to quantities in 
e,PQ~p' ~ x^)' shift to a coordinate u such that u = e^~''<3, m^, = e^'^f^. Defining Xq 
as the 't Hooft coupling fl2.39p at the quark mass scale, inserting the value of Qc in (12.11) 
and defining Tq = {e'^^'^y/—gtt9xx) \p=pq ^-s the tension of a hypothetical fundamental string 
stretched at constant p = pg, we can write the estimate for the meson masses as: 




n 




(2.60) 




(2.61) 
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It is important to stress once again that this expression is written only in terms of IR 
quantities, once we discard terms of order = e^''"^*, namely contributions like logu* 
have cancelled out. Notice that the upper limit of the integrals can be taken to infinity if 
we again insist in discarding 0{u~^) contributions. The expression ( 12.6ip is a neat example 
of how, even having a Landau pole, the holographic set-up is able to consistently obtain IR 
predictions, in exactly the same spirit as in field theory. We can perform numerically the 
integration in fl2.6ip . and we get [70] : 

1 

M^") ^ (5.2 - 6 X 10^=^^-^ + . . . ) , (2.62) 

where in order to substitute eg as in fl2.40p we have used Vol{X^) = ^vr^, Vol{X^) = ^ir^. 
The expression (12.62^ is the result quoted in ^70j, apart from a different factor of 2 in the 
definition of Xq. 

The lesson we want to take from this section is that there is a well defined method to obtain 
the shift produced by the flavor quantum effects on the meson spectrum (or, eventually, on 
any physical observable) as an expansion in the parameter e ~ A Nf / which weighs the 
flavor loops. Heuristically, it may be useful to think of the computation leading to (12.62^ as 
(partially) a strong coupling analogue of the Lamb shift corrections of QED. 



2.5 Black hole solutions: D3-D7 quark-gluon plasmas 

In this subsection we will review the results in [68]. We start by showing how one can find 
a black hole solution which includes the backreaction effects due to massless quarks. To 
perform this analysis it is more convenient to work with a new radial variable a such that 
the metric takes the form: 

ds^ = h-^ [-b de + dxl] +h^[h e^»+2^ da^ + e"' ds^j, + e'^ {dr + Ake? ] ■ (2.63) 

Notice that we have introduced a new function h which parameterizes the breaking of Lorentz 
invariance induced by the non-zero temperature T. All functions appearing in the metric 
fl2.63p . as well as the dilaton 0, depend on a. Moreover, the RR field strengths F5 and Fi 
are given by the ansatz fl2.33p with the function p = 1. We remind the reader that fixing 
p = 1 corresponds to taking massless quark^. 

In this non-supersymmetric case we will not have the first-order BPS equations at our 
disposal and we will have to deal directly with the second-order equations of motion. Ac- 
tually, since all the functions we need to compute depend only on the radial coordinate a, 
it is possible to describe the system in terms of a one-dimensional effective action. One can 

^In reference [68' , the more involved case of massive quarks p 7^ 1 was also discussed. An extra compli- 
cation is the necessity of finding the non-trivial D7-brane embeddings in the backreacted geometry. 
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find this effective action by directly substituting the ansatz in the gravity plus branes action 
ffLHD . One gets: 

+ 24 6e2^+^5-46e^^+^5 - - ^he^^+^^ - 4Q^ 6 e^^+^^+^z^ ^ ^2.64) 

In fl2.64p Vi^3 denotes the (infinite) integral over the Minkowski coordinates. The second 
derivatives coming from the Ricci scalar have been integrated by parts and, as is customary, 
only the angular part of is inserted in the F| term (otherwise the Qc would not enter the 
effective action since, on-shell, = due to the self-duality condition). The last term in 
f l2.64p . proportional to Q/, comes from the DBI contribution in (11.41) . Notice also that the 
WZ term does not enter (I2.64p because it does not depend on the metric or the dilaton. 
The equations of motion stemming from the effective action fl2.64p are: 

9^ (log 6) 
djilogh) 

dig 
dlf 



dl(j) = gjfoe^'^+ss +4Q^5e<^+6<;+2/ . (2.65) 

It is straightforward to check that these equations solve the full set of Einstein equations 
provided the following "zero-energy" constraint is also satisfied: 

+ ^-^[^ + 8d^g + ^d^f) - 24 6e2^+^^ + 46e^^+^^ + 

+ 91^ + 9lbe''^+'^ + 4Qfbe^+^^+^f . (2.66) 

This constraint can be thought of as the aa component of the Einstein equations or, 
alternatively, as the Gauss law from the gauge fixing of g^a in the ansatz (12.631) . 
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Let us now find a solution of the system of equations fl2.65p and of the "zero-energy" 
constraint f l2.66p that corresponds to a black hole for the backreacted D3-D7 system. We 
will require that such a solution is regular at the horizon and tends to the supersymmetric 
one at energy scales much higher than the black hole temperature T. Actually, the biggest 
advantage of the radial variable a introduced above is that the equations of motion of h and 
h in (I2.65P are decoupled from the ones corresponding to the other functions of the ansatz. 
These decoupled equations can be easily integrated in terms of an integration constant r/^ as 
follows: 



6 = e^'^>, h = ^A\-e<"). (2.67) 



where a G (— C)0, 0). We now define a new radial coordinate r by means of the relation: 

r 

Then, b and h take the form: 



e4-> = l-!A, re{rf„+oo). (2.68) 



1 - ^ , h = (2.69) 



with = QJA. Notice that h is given by the same expression as in fl2.50p . Moreover, it is 
clear from (12.690 that r = r^is the position of the horizon and, thus, the extremal limit is 
attained by sending to zero. In terms of r the metric takes the form: 



it^ L V r^/ J r"^ 1 _ 

J- ^4 

(2.70) 

where we have defined the functions / and g as follows: 

pf p9 

= — , = - . (2.71) 

r r 

In order to determine completely the background we still have to solve eqs. (I2.65P and fl2.66p 
for /, g and the dilaton (f). We will find this solution by introducing a reference UV scale r* 
and by expanding the functions in terms of the parameter e^, defined in fl2.38p . We will impose 
that the functions /, g and (j) are equal to the SUSY ones of fl2.52p when the extremality 
parameter Vh vanishes. Moreover, we will also require that these functions coincide with 
those in (I2.52p at the UV scale r^. These conditions fix uniquely a solution of (12.650 and 
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fl2.66p . Up to second order in e* this solution is given by: 




(2.72) 



where Li2{u) = J2'^=i ^ is a polylogarithmic function. The functions written in fl2.72p 
determine a geometry that is regular at the horizon r = r^. In the next subsection we will 
study its thermodynamics and we will extract some consequences for the dual field theory 
with dynamical quarks at non-zero temperature. 

Let us conclude this section with some comments on the stability of our perturbative 
non-extremal solutions. A possible way to check for the latter is to consider worldvolume 
fluctuations of a D7-brane in the setup. If, as in our cases, the brane corresponds to massless 
flavors, the related quasi-normal modes on the unflavored background all have frequencies 
with a negative imaginary part of the order of the temperature, signaling stability [75], [76] . 
This result cannot be changed in the flavored case when a perturbative expansion in is 
done. Thus, in our regime of approximations, stability with respect to those fluctuations is 
guaranteed. 

2.5.1 Thermodynamics of the solution 

In the previous subsections we have defined the backreacted background in terms of an 
arbitrary UV scale r^, as an expansion in powers of the parameter e^, written in (12.401) . This 
scale r^, should be well separated from the Landau pole scale in order to avoid having the 
pathologies of the latter. Moreover, we are interested in analyzing the physical consequences 
of this background at energies much lower than the UV scale r^,. In a black hole background 
dual to a quark gluon plasma the natural IR scale is the location of the horizon, which 
should be related to temperature T of the plasma. Accordingly, we define eu as: 




(2.73) 
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where, in what follows, the subscript h means that the quantities are evaluated at the horizon 
r = Th- Thus, \h is naturally identified with the 't Hooft coupling at the scale of the plasma 
temperature. We therefore have: 



e. + e:iog^ + 0(e^) 



(2.74) 



We will use this relation to recast the expansions in powers of e* as series in eh- We will 
assume in what follows that r/j is well below the reference UV scale r* to ensure that the IR 
physics does not depend on the UV completion of the theory. In a Wilsonian sense of the 
renormalization group fiow, the UV details of the theory should not affect the IR physics. 
Moreover, since, as we will see below, Vh is proportional to the temperature (at leading 
order), we have: 

^ = ^ + 0(eD, (2.75) 

and T[d\h/dT) = ehXh at leading order. These relations refiect the running of the gauge 
coupling induced by the dynamical fiavors. 

The thermodynamic properties of the black hole solution are determined by the metric 
functions at the horizon. After neglecting terms suppressed in powers of r^/r^, the values of 
the functions / and g at r = Vh can be obtained from (12. 72 p . One gets: 



ofh 



17 



24 ^ 1152 



o9h 



1 + 



24 128 



+ 0(e 



(2.76) 



The black hole temperature can be obtained by requiring regularity of the euclideanized 
metric and by identifying the temperature with the inverse of the period of the euclideanized 
time. A simple computation yields: 



T 



2rh 



1 - 5- -si^' 



(2.77) 



where in the last step we have used the values of fh and gh written in (12. 76 p . 

The entropy density s is proportional to A^, the volume at the horizon of the eight dimen- 
sional part of the space orthogonal to the t, r plane (where i is the Euclidean time), divided 
by the infinite constant volume of the 3d space directions V3. From the general form of the 
metric we get that: 



2n As _ rlR^e^^h+k Vol{X^ 



TT 



2Vol{X^ 



-n: 



TT 



1 19 



which in terms of the temperature reads: 



IT 



2Vol{X'^ 



(2.78) 
(2.79) 
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As for the other thermodynamic quantities which will follow, the leading term of this 
formula is the well-known unflavored result. The 0{eh) term was already calculated in [76] 
with the probe brane technique, in the = case. Here we have re-obtained this result in 
a quite standard way, by computing the increase of the horizon area produced by the flavor 
branes. This can be considered as a crosscheck of the validity of the whole construction. 
Finally, the order was first obtained in [68] . 

The ADM energy of the solution can be computed as an integral of the extrinsic curva- 
ture of the eight-dimensional hypersurface of constant time and radius. This calculation is 
straightforward and has been done in the appendix B of [68j, with the result: 



Eadm 



Vs 8Vol{X^ 



1 + ^e,iT) 



1 



e,(T)2 + 0(e,(T)= 



(2.80) 



Again, terms suppressed as powers of ^ have been neglected. Moreover, since in the following 
derivatives with respect to T are going to be taken, we find it convenient to make explicit 
that eh depends on T (see fl2.75p ). Eq. f l2.80p yields the energy density of the plasma and, 
thus, it allows us to study the full thermodynamics. Indeed, from f l2.75p and f l2.80p we get 
immediately the heat capacity (density): 



- — N^T^ 

2Vol{X^) " 



l + ^6,(T) + lie,(T)2 + 0(e,(T)= 



The free energy density, and so (minus) the pressure, reads: 



F 

"^3 



T7 = -P 



Ts 



TT 



8 Vol{X'- 



1 + -^^KiT) + -eh{Tf + 0{eh{Tf) 



(2.81) 



(2.82) 



Notice that, consistently, this satisfies the relation s = dxp (where it is crucial to take f l2.75p 
into account). This result is confirmed by the direct computation of F from the renormalized 
euclidean action (see, again, the appendix B of [68]) and also by the calculation in [TT] of 
the correlator of the tensorial mode in the hydrodynamical approximation. 
The speed of sound Vg is obtained by combining f l2.79p and fl2.8ip . namely: 



s 

Cv 



(2.83) 



Note that the correction to the speed of sound, which measures the deviation from confor- 
mality, only appears at second order and that the sign of the correction is consistent with 
the bound vl <^ conjectured in [7^. It is also interesting to point out that the solution pro- 
vides a direct measure of the breaking of conformality at second order, namely the so-called 
interaction measure, given by: 



£ — 3p 

^ lQVol{X^ 



(2.84) 
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Let us now analyze the viscosity of the plasma predicted by the flavored black hole. Since 
we are not introducing higher derivatives of the metric in the action, the usual theorems apply 
and the shear viscosity rj saturates the Kovtun-Son-Starinets bound [78], i.e. rj/s = l/An. 
Therefore, the shear viscosity 77 can be obtained by dividing by An the entropy density 
written in f l2.79p . Again, the first-order term coincides with the one calculated in the probe 
approach in [79] while the second-order result was first computed in [68] . On the other hand, 
one can also compute the bulk viscosity ( for this model, with the result [7T] : 

C = 72Vd{X^) ^'^^' + 0{e,{Tf)] . (2.85) 

Interestingly, the value of ( written in f l2.85p saturates the bound proposed in [80] : 




(2.86) 



For the computation of other transport coefficients, we refer the reader to [71] . 



2.5.2 Energy loss of partons 



One of the main phenomenological applications of holography is the analysis of the energy 
loss of a parton that moves through a quark-gluon plasma. One of the measures of this energy 
degradation is the so-called jet quenching parameter q, which is a transport coefficient that 
measures the bremsstrahlung experienced by a parton probe due to its interactions with 
the quarks and gluons of the plasma [Hlj. At very high energy, and using the eikonal 
approximation, the authors of [82] found a non-perturbative prescription for calculating q as 
the coefficient of in an almost light-like Wilson loop with dimensions ^ L. By using 
the generic formula in [83] (and cutting the integral at r^,), we can write: 



71 a 



'9r 



Qxx a/ 9 XX + Qtt 



zdr 
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(<t'-4>h) 
2 



dr 



(2.87) 



The dilaton enters the formula because we are considering the Einstein frame metric. By 
plugging in f l2.87p the expressions of f, g and written in (12.72 p . and by performing the 
corresponding integrals in r, one gets g as a power series expansion in e*. In the course of 
this calculation we will neglect terms that are suppressed by powers of r/^/r* and we will 
write the result in a series in eh rather than in e^,. In terms of gauge theory quantities one 
gets [68]: 



TT 



rp3 



1+ (2 + n)eh + iel + 0{el) 



where we have introduced a constant 7: 
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96 ^ 48 
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-C + — 4F3 
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0.5565 



(2.88) 



(2.89) 
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with C ~ 0.91597 being the Catalan constant. Notice that the flavor correction to q is 
positive, i.e. fundamentals enhance the jet quenching. Actually, eq. f l2.88p can be used to 
estimate this enhancement in the extrapolation to the realistic RHIC regime. Let us take 
= 5^ N^ = Nf = 3 and a, = g^M/^^ ~ 1/2- Then, A,, ~ Gtt and e;, ~ J ~ 0.24. 
Using this value in f l2.88p we have that q is increased by 20%. For example, at T = 300 
MeV we obtain q ~ 5.3 (Gev)^/fm, to be compared with the value [82] q ~ 4.5 (Gev)^/fm 
of the unflavored plasma (the RHIC values are g ~ 5 — 15 (Gev)^/fm). It is also interesting 
to rewrite f l2.88p in terms of the entropy density s. One gets: 



which shows a deviation (driven by eh) from the general expression put forward in [81]. In 
this setting, the presence of fundamentals and the breaking of conformality are inevitably 
mingled. It would be interesting to have the dual of a conformal theory with fundamentals 
to check whether the conjecture [H] holds in such situation. That was analyzed in [57] in 
a non-critical string framework and, interestingly, the result differs from \8^. The caveat is 
that the model studied in [57] suffers from the usual problem of gravity-like approaches to 
non-critical strings, namely there are uncontrolled approximations. 

Another way of characterizing the energy loss of a parton probe in the plasma is by 
modeling it as a macroscopic string attached to a probe flavor brane. The string is dragged 
by a constant force / which keeps its velocity v flxed and transfers to the parton energy and 
momentum, which is then lost in the plasma at a constant rate. This energy loss is measured 
by the drag coefficient /i, which relates the force / and the parton momentum p: f = fip. To 
compute this drag force one can apply the general procedure of refs. [851 ESI [87]. By using 
the Nambu-Goto action for a string in the black hole background one gets that the rate of 
momentum transferred to the medium is given by 



dp 1 ^ rl -fire) V V 

-T- = C = , = -a Mkin , , 2.91 

dt 2Txa' 2Tia' ^JY^ ^JY^ 

where C is the constant determined from the equation gxx{j'c)9tt{j'^ -fC^ = with the point 
Tc given by guiy,!) + gxx{i^c)v'^ = 0, namely Vc = rh{l — f^)"*. In f l2.9ip we have introduced, 
following ref. [85], the kinematical mass M^in such that p = Mjsi„-^=^. From fl2.9ip . using 
(EUD, (frZlD, (ETTD, we find: 



a Muin = : I 

^ 2 ^Vol{X^) 



l + ^(2-log(l-t;2))e,+ (2.92) 



1 



[44 - 20 log(l -v^) + 9 log^(l - f 2) + 12Li2{v^)] el + 0{e 



3^ 

As happens with the jet quenching, the energy loss (at flxed v) is enhanced by the presence of 
fundamental matter. The quantity /i M^j^ grows when increasing the velocity. From 
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formally, it would diverge as f — )■ 1. However, (12.921) is not applicable in that limit since we 
have to require eh log(l — f ^) ^ 1 for the expansions to be valid. 

2.6 A discussion on the range of validity 

We now discuss, following ^68j, the restriction on the physical parameters needed for the 
deconfined flavored plasma solution to be physically meaningful. Before we go on, two 
comments are in order: first, notice that, even if we will use here the plasma temperature 
as the IR scale at which the relevant physics takes place, this can be substituted by any 
other IR scale, depending on what one wants to study. Thus, for instance, when computing 
meson masses at zero temperature as in section l27i| the discussion below holds, just taking 
into account the IR scale there is set by the quark mass. Second, notice that the restriction 
of small e=K (which leads to Nf <^ N^) of the D3-D7 case at hand comes from the existence 
of a Landau pole. In holographic theories in which there is no Landau pole in the geometry 
(sections El H]), there is in principle no restriction to Nj. In particular, it is possible to 
consider in those theories Nf to be of the same order as Nc. 

As we have already remarked, having a pathological UV means that there must exist a 
separation of scales between IR and UV. Concretely, there must exist a hierarchy, which in 
terms of the r radial coordinate reads: 

r/, < < r„ < r^p . (2.93) 

The quantity Vh sets the scale of the plasma temperature ~ Kjr ~ T, which is the 
scale at which we want to analyze the physics. The point r^p is where the dilaton diverges, 
signaling a Landau pole in the dual theory. At a scale the string solution starts presenting 
subtler pathologies, whose discussion we delay until the end of this section. Finally, sets 
an (arbitrary) UV cutoff scale ~ ^uv- The solution (I2.72p will only be used for r < r*. 
In a Wilsonian sense of a renormalization group flow, the UV details should not affect the 
IR physical predictions. This feature is reflected in the fact that physical quantities do not 
depend (up to suppressed contributions) on r^, or functions evaluated at that point, but only 
on IR parameters. Even if the precise value of is arbitrary, we have to make sure that it 
is possible to choose it such that it is well above the IR scale (so that the UV completion 
only has negligible effects on the IR physics) and well below the pathological r^, r^p scales 
(so that the solution we use is meaningful and the expansions do not break down). To this 
we turn now. 

Let us start by computing the hierarchy between r* and tlp. Since at r^, we can approxi- 
mate the solution by the supersymmetric one, we can read the position of the Landau pole 
from fl2.45p . If we insert the approximate relation between radial coordinates r f» a/o'c'', we 
find: 

— ^ < 1 , (2.94) 

rLP 
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as long as e* ^ 1. 

Moreover, one has to make sure that the Taylor expansions 
Vfi < r < r^:. This of course requires e,,, ^ 1, but also that e^, 



I2.72P are valid in the region 
log — ^ 1 (notice that the 



absolute value of the logarithm can be big because rh ^ r^:). This means that ^ ^ e~~ . 
On the other hand, when computing physical quantities in the previous sections, we always 
neglect quantities suppressed as powers of ^ ~ -j^^. This is the order of magnitude of the 
corrections due to the eventual UV completion of the theory at r*. One has to make sure 
that the corrections in e^, we are keeping are much larger than the neglected ones, namely 
e=K ^ ^. In summary, we have the following hierarchy of parameters (in the following, in 
order to avoid overly messy expressions, we insert the value of e^, for the = case, 
remembering that for a generic X^, its value is given by ( I2.40p ): 



1 



7r2 JV^ 



e ^* ~ e ^* ^/ < — ~ < ~ -^f^ < 1 . (2.95) 

As long as e^, ~ ^^r^' ^ ^' there always exists a range of r^, such that this inequality 

is satisfied. Since we focus on the IR physics of the plasmas, at the scale set by their 

temperature, the actual physical constraint on the parameters will be <^ 1 , which we 

have written in terms of the coupling at the scale of the horizon, Xh = A*(l + 0{e^:)). 

On top of this, we have to make sure that the SUGRA+DBI+WZ action we are using 

is valid. As usual, the suppression of closed string loops requires A^c ^ 1 whereas the 

suppression of a'-corrections is guaranteed hj Xh ^ 1. We have written the D7-brane 

worldvolume contribution to the action as a sum of Nf single brane contributions. This is 

justified if the typical energy of a string connecting two different branes is large (in a' units). 

1 

Since the branes are distributed on a space whose size is controlled by _R ~ a/o', we again 
need A/j ^ 1. The smearing approximation will be good if the distribution of D7-branes on 
the transverse space is dense, i.e. Nf ^ 1. The discussion up to now is summarized in the 
following validity regime: 

AT, > 1 , A^ > 1 , AT; > 1 , eh = < 1 . (2.96) 

Finally, we want to find the regime of parameters in which the flavor corrections are not only 
valid but are also the leading ones. With this aim, we ought to demand that the leading 

_ 3 

a'-corrections to the supergravity action (which typically scale as A^ ^ due to terms of the 
type a'^TZ"^) are smaller than the flavor ones, controlled by eh, namely: 



3 



Xh~'<^eh. (2.97) 

Demanding that corrections to the D7-branes contributions (for instance curvature correc- 
tions to the worldvolume action itself or corrections produced by possible modifications of 
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the brane embeddings due to curvature corrections to the background metric) are subleading 
does not impose any further restriction. The reason is that their contribution is typically of 
order ehXj'^'^ for some c > which is always subleading with respect to et as long as fl2.96p is 
satisfied. 

The holographic a-function 

As discussed in [TOj and mentioned above, the string solution starts presenting pathologies 
at a scale < r^p, where the holographic a-function is singular. The utility of the solution 
for r > is doubtful, but since we have only used the solutions up to r^: <^ in order 
to derive the IR physics, this subtlety does not affect the physical results. We now briefly 
review the argument in [70], which used the backreacted Klebanov-Witten solution at zero 
temperature. The qualitative picture holds for the rest of the cases addressed in the present 
section and for the case of section |5] too. 

Let us start by considering the metric of a generic dimensional reduction to five dimensions, 
giving a 5d Einstein frame metric of the form (the u here is, obviously, a redefined holographic 
coordinate, namely u = u{r)): 

dsl = H{uf'^ [rfx? 3 + P{u)du^] . (2.98) 

In standard set-ups, the function H(uy^^, which can be roughly identified with the dual field 
theory energy scale, monotonically varies with the radial coordinate. This is also required 
in order for the "holographic a-function" [88] : 

a{u) ~ /3{uf/^H{uy/^[H'{u)]-^ , (2.99) 

to be finite]^ Instead, the function H{u) is not monotonic here: it increases with u from zero 
up to a maximum at a point Ua and then it decreases back to zero where h vanished. In 
the flavored supersymmetric KW case, the H and a functions read: 

H{p) ~ he'f+'^ , a(p) ~ h'/' e'f H'/'[d,H]-\ (2.100) 

where we have not written unimportant overall factors. A representative plot is given in figure 
|5l The non-monotonic behavior of H implies that the holographic a-function is singular and 
discontinuous at the "a-scale". From the plot, we see that ~ e^" is below, but not 
parametrically separated from rip. 

®The monotonicity of H{u) also plays a crucial role in holographic computations of the entanglement 
entropy, see [89]. The notations of that paper are used in the equations above. 

For the present discussion and in particular for figure [S] we will choose the additive integration constant 
of h such that h is zero at the Landau pole. The specific point Ua (namely Tq) at which this UV pathology 
sets in depends on this choice. Again, we stress that the important point is the IR results do not depend on 
this choice (modulo suppressed contributions) as long as 3> . What we show here is that the integration 
constant can be naturally chosen such that this condition is satisfied. 
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Figure 5: The function H in the massless-flavored KW model at zero temperature. 

3 A dual to A/" = 1 SQCD-like theories 

In the following section, we will study a system that in some sense is qualitatively different 
from those of the previous sections, though the procedure to deal with the addition of flavors 
is identical. The main qualitative difference will be that there need not be a hierarchical 
difference between the number of flavors and the number of colors. The case treated here will 
represent the addition of fundamental matter to a field theory that is originally confining and 
four dimensional at low energies, but that gets some higher dimensional completion in the 
UV (in principle this allows one to extend the range of the radial coordinate to arbitrarily 
large values). Some of the qualitative changes that observables of a confining theory undergo 
when fundamentals are added will be discussed. The developments described in the present 
section were applied to model possible aspects that could appear in physics beyond the 
standard model, as we will briefly mention below. 

More concretely, in this section, we will study a dual to a version of A/" = 1 SQCD. The 
model is based on D5-branes wrapped on two-cycles inside the resolved conifold-leading to 
a geometry related to the deformed conifold. We will first briefly present the model without 
flavors, then study the addition of flavors following the ideas described in the first sections of 
this paper: kappa symmetric embeddings, smearing, backreaction, system of BPS equations, 
particular solutions to this system and finally present a set of checks that the correspondence 
we are proposing is valid and robust, we also explain some predictions about the field theory 
obtained with the string background. 

3.1 The model without flavors 

The proposal is to construct a dual to a field theory with minimal SUSY in four dimensions 
using wrapped branes. Ideas of this kind were first explored by Witten in the early days of 
AdS/CFT. In the paper [90], Witten presented a model dual to a version of Yang-Mills theory 
(with an extra massless scalar and UV-completed by an infinite tower of massive vectors, 
scalars and fermions), by wrapping a set of Nc D4 branes on a circle with SUSY-breaking 
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boundary conditions. 

The idea here is very similar, only that we will work with D5-branes and we will preserve 
some amount of SUSY. We will compactify the five branes in a very subtle way (involving a 
twisting of the 6-d theory) so that only four supercharges will be preserved in the compactified 
theory for all energies |9lj"1 . in other words, the partial SUSY breaking is not due to the 
presence of relevant operators, like mass terms. This kind of compactification of the six 
dimensional theory living on a stack of D5-branes (when the D5's wrap a two-cycle inside 
the resolved conifold) was well studied in various papers, see [93] for various reviews. We 
will follow mostly the detailed study of [94] . 

One can show that a very generic string background describing a stack of Nc D5-branes 
wrapping a two cycle and preserving four supercharges includes a metric, RR-three form 
-F3 = dC2, a dilaton (f){p) and is given by. 
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((wi + a{p)d6y + {002 — a{p) smOdipY^ H — {003 + cosOdipY 



(wi + h{p)d6) A {u2 - h{p) ^wiBdip) A (^3 + cos 9d(f) + 



b'dp A {-d6 Aui + sin 9dip A U2) + (1 - Kp^) sin 9d9 A dip A 



where Ui are the left-invariant forms of SU(2) 

Ui = cos tljd9 + sin sin 
U2 = — sin tpd9 + cos ip sin 
U3 = dip + cos 9dip>. 



(3.1) 



(3.2) 



For convenience, below we will set the parameters Qs = a' = 1. The presence of the Nc color 
D5-branes is indicated in F3 that satisfies the quantization condition: 



2 '^(10) Js^ 
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NrT, 



c-l-5 ■ 



(3.3) 



The S"^ on which we integrate is parameterized by 9,(f,ip. 

We then impose that a fraction of SUSY is preserved, hence we need to impose some 
projections on the Type IIB spinors and a set of BPS equations refiecting this arise (see the 
appendix B in the paper [11] for generous details). The BPS equations are non-linear, first 



"'^"'^The fact that a twisting procedure (see [92] for a very nice presentation of this idea) is at work, imphes 
that even in the far UV, the theory is still preserving only four supercharges . 
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order and coupled for the functions of the background in eq. fl3.1l) -for details see appendix 
B in [llj. Certainly, solving first order equations is simpler than solving the second order 
Einstein equations; nevertheless the BPS equations for the functions (0, h, g, k, a, b) are non- 
linear and coupled, rendering the problem complicated. 

It is technically convenient to make a 'change of basis' to another set of functions, so that 
the BPS equations become first order and nonlinear (of course) but can be decoupled, and 
then solved independently. A change of variables that does the job partially was obtained 
in [21]. The change of basis is from the set of functions [(j),h,g,k,a,b] into the functions 
[P, g, r, $, Y, a]. The map reads [24], 



2^ If P^-Q^ \ 20 „ V- ^ 2k Psinhr , a 

e = 7 7^ , ^ , = Pcoshr - Q, e^" = AY, a 



and 



4 Vi^coshr-gy ' ' Pcoshr -g' A^, 

(3.4) 

$ = (p2_g2)yye2^. (3.5) 

As explained in detail in [23] (see section 3 of that paper), the BPS equations can be solved 
one by one for these new functions, obtaining: 



g(p) = (g„ + Nc) coshr + A'c(2pcoshr - 1) 
1 

sinh(2p — 2po) 



sinh r(p) = , -, coshr(p) = coth(2p — 2po) 



P' 

Yip) = IT' 



e 



4^ _ e^^°cosh(2po) 



(p2-g2)Fsinh'r 



a = tanhT{Q + N^) = — —77^ TT^- 3.6 

smh(2p — 2po) 

Note that both Q and the dilaton are given algebraically in terms of the rest of the functions 
parametrizing the backgrounds. Here Qo and (po are constants of integration and we have 
chosen the integration constant in the dilaton field equation 0o such that it admits a smooth 
limit as po — )• —00 (this limit gives t = a = and so corresponds to what in the paper [95] 
were called type A backgrounds). 

The function P satisfies the following second order equation: 

We will refer to this equation as the 'master' equation, since once we have a solution of (13. 7p 
all other functions are determined via (13. 6p . 
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3.2 Some solutions 



There are many solutions to the master equation (13. 7p . A very simple one is given by 

P = 2N,p, Qo = -N,. (3.8) 

Once processed back, one computes the functions in the original background of eq. (13. ip 
and one recovers an old solution [96]. To avoid nasty singular behaviors, in the following, 
we will choose the value of the integration constant Qo = —N^, so that the first term in the 
expression for Q{p), namely {Qo + Nc), vanished. 

Aside from the simple solution presented above, there are a variety of very interesting 
solutions. For example when the function P{p) is given near p = by the following Taylor 
series: 

where hi is again an arbitrary constant (notice that for hi = 2Nc we get back to the solution 
in eq. (13. 8p : we will also assume that hi > 2Nc). It is interesting that this solution can 
be numerically connected in a smooth way with a solution for large values of the radial 
coordinate (p — oo) that differs greatly from the linear behavior of the solution in eq. (13. Sp . 
In this case, it is given by 



P ~ et^ 



c(l - |pe-^0 + 7:^(256p2 + 256Qop + UANj + QAQDe-^ + ©(e-^^l . (3.10) 
3 d4c J 



These solutions were studied explicitly in section 8 of the paper [TT] and have a variety of 
interesting applications that we will briefly mention in the following sections. 



3.2.1 An exact recursive solution 

There is one recursive way of obtaining solutions, described in [97], that basically uses the 
fact that the master equation ( 13. 7p can be written as (we choose here and in the following 
Po = 0) 

dp (s{P^ - Q^)P') + AsP'QQ' = 0, s{p) = sinh^ r = (3.11) 

smh (2p) 

Integrating eq. (I3.1ip twice we obtain 

P^ - SQ'^P + Q r dp'QQ'P + 12 r dps-^ T dpsP'QQ' = c^R{pf, (3.12) 

J P2 J p2 J pi 

^^If we do not make this choice, the space ends before p ^ po, since Q > P possibly giving place to 
geodesically incomplete spaces and a divergent dilaton. Hence, we will choose the term proportional to 
coshT(p) in Q{p) to vanish. 
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where 

R{p) = (cos^a + sin3a(sinh(4p) - Ap)Y'^ , (3.13) 

being (c, a) the two integration constants of the master equation. 

Following [2j] we write P in a formal expansion in inverse powers of c - the integration 
constant encountered above - as: 

oo 

P = Y,c'^"Pi-n. (3.14) 

n=0 

Inserting this expansion in eq. f l3.12p we obtain recursively 

Pi = R, 

Po = 0, 

P-1 = -ipr' ( -3Q'Pi + 6 r dp'QQ'P^ + 12 r dp's-^ r dp"sQQ'P[ 



n+2 / n— m+3 



P-n-2 = ~'^Pl ^ i ( '^P^P^-'rn.Pm-n-2 + Pl-mPl-kPm+k-n-2 1 " "iQ^ P-n 

lm=l \ fc=l / 

+6 r dp'QQ'P.n + 12 r dp's-^ T dp'sQQ'PL^ I , n > 1. (3.15) 

Jp2 "'P2 ■'Pi J 

It follows by induction that = for all even k. The large p expansion of these solutions 
coincide with that described in eq. fl3.10p Once again, solutions written in this form have 
interesting applications to the physics of cascading quivers on the baryonic branch pS] , 
[nS]. We will not study the physics encoded in the solutions described above, suggesting the 
interested reader to consult the papers [93], [11], [97] . 

There is another set of solutions, proposed in [57] and whose physics content was developed 
in [100j . [T0T] that correspond to what are called "walking solutions". The idea here is to 
construct string backgrounds such that the dual QFT has a gauge coupling with very slow 
running (or "walking" coupling). See the papers [^. [lUUj . llUl] for detailed explanations on 
the physical implications of these solutions. 



3.2.2 A comment about the dual field theory 

The system of wrapped D5-branes has a field theory realized on its worldvolume, whose 
dual background and various solutions were described above. The field theory is a version 
of minimally SUSY Yang- Mills. Again, some UV completion takes over at high energies 

""^^We are not saying that Super- Yang-Mills needs a UV completion, just that the system of D5-branes 
realizes a theory with these characteristics. 
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The field theory is minimally SUSY (A/" = 1) and its perturbative spectrum, aside from a 
massless vector multiplet contains a tower of massive vector and chiral multiplets. A careful 
study of the perturbative dual field theory obtained by compactification and twisting of the 
six-dimensional theory living on (unwrapped) D5-branes was done in [94]. In that paper, the 
degeneracies and masses of the (perturbative) states in the tower mentioned above are given. 
More interestingly, the authors of [91] showed that the theory is equivalent to J\f = 1* Yang- 
Mills in a particular Higgs vacuum, where the extra dimensions appear by deconstruction. 
In this sense, we will think of the theory without flavors either as a six-dimensional theory 
compactified or as a four-dimensional theory with an infinite set of fields. 

For our purposes, it will be enough to use the fact that the lagrangian of the field theory 
reads 

L = Tt[~FI - tX^D.X + Wk, W)] , (3.16) 

where and Wk represent the infinite number of massive chiral and vector multiplets and 
W denotes the massless vector multiplet. The term L($fc, Wk, W), represents all the kinetic 
terms and interactions that can be deduced from [9l] . More comments about this field theory 
can be found in the appendix A of p5j. 

In what follows, we will summarize the procedure of adding flavors to this field theory. 
The flavor branes, in this particular case are D5-branes. 

3.3 Addition of flavors 

The study of supersymmetric embeddings in backgrounds of the form of eq. (13. ip . more 
precisely for the solution given in eq. (13. 8p was initiated in the paper [18]. There the 
eigenspinors of the kappa symmetry matrix were found to be the spinors preserved by the 
background for a variety of D5-brane embeddings. For the purposes of this review, we 
will focus on the "cylinder embeddings" described in section 6.3 of [18] and in more detail 
in section 6.5.3 of the third paper in reference [93]. In this case the flavor D5-branes are 
extended along the R^'^ Minkowski directions, on the radial direction p and also wrap the R- 
symmetry direction ip. Intuitively, the flavor branes are localized in the directions {6, 6, (f, (p), 
but interestingly enough, any constant value of these coordinates ensures that we have a 
kappa symmetric configuration. This is a very important fact, as we can put one flavor 
brane "at each point" of the four manifold S[^, 9, (p, cp] and still have a SUSY configuration. 

This is precisely what we will take advantage of when smearing. Let us see this in more 
detail: if, as discussed in the first section, we write the action describing the closed strings 
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(IIB) and the open strings (BIWZ), we will have 



2 '^^10) 



S = I d^\^ 



Nf Nf 

-T5 V / rfSet,/^ + ^5 V / P[C,] , (3.17) 

where the integrals are taken over the six- dimensional worldvolume of the flavor branes Aie, 
and ^(6) stands for the determinant of the pull-back of the metric in such a worldvolume. 

As discussed in previous sections, we then think of the Nf — )■ 00 branes as being ho- 
mogeneously smeared along the four transverse directions parameterized by the coordinates 
6, (fi and 6, (p. The smearing erases the dependence on the angular coordinates and makes it 
possible to consider an ansatz with functions only depending on r, enormously simplifying 
computations. One has: 

Nf 

-^sX]/^ rf^xetyi^ J rfi^sin^sin^^ety^, (3.18) 

Nf 

r 1^ N r 

^[^6] ^ y yoKyA)^c^,) , (3.19) 

where we have defined Vol{y^ = sin 6 sin 6d6 A dip A d6 A dip and the new integrals span the 
full space-time. We will need the expressions (with the choice explained above a' = Qs = I). 

^5 = 77^ , 24o. = (27r)^. (3.20) 



From here, we will have a set of BPS equations describing the dynamics of this open-closed 
string system. The same change of basis with the purposes described around eq. (13. 4 p can 
be performed - see [M] for details. The solution in this case is dependent on the number of 
flavor branes Nf and reads (reinstating momentarily the integration constant Po): 

sinhr = . , / rr, (3.21) 

smh(2(p - po)) 

for the function r, while for Q, $ we have, 

^ 2Nc-Nf\ , 2N^-Nf, , , , , 

Q= iQo + — — ^jcoshr + ^^^ ^(2pcoshr-l), (3.22) 



^A{^-4>o) _ cosh2(2p. 



(3.23) 
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In the case with flavors, hke in the unflavored case previously discussed, both Q and the 
dilaton are given algebraically in terms of the rest of the functions parametrizing the back- 
grounds. As before, po, Qo and (po are constants of integration and we have chosen the 
integration constant in fl3.23p such that it admits a smooth limit as po — )■ — oo (this limit 
gives r = cr = and so corresponds to the type A backgrounds). The function Y is deter- 
mined in terms of P as 

Y=^-{P' + Nf), (3.24) 

while the only remaining unknown, the function P, then satisfies the new decoupled second 
order master equation 

P" + iP' + N,) {^^-^^^ + '^^^^ - - 2p„)) ^ 0. ,3.25) 

One can redefine P{p) = Ncp{p) and factor out Nc from the master equation. We will 
mention some solutions to eq. f l3.25p . that explicitly include the quotient x = hence the 
solutions will capture the nontrivial physics of the fields transforming in the fundamental 
representation of the gauge group. 



3.4 Study of solutions 

We now describe various solutions to the "flavored" master equation (13.251) . Some solutions 
were found exactly, for the particular relation Nf = 2Nc while some other are known as 
asymptotic expansions, near the UV (large p) and the IR (small p). In these latter cases, a 
smooth numerical interpolation can be found. 



3.4.1 Exact solutions for Nf = 2N^ 

One can find some exact solutions for the case Nf = 2Nc or x = 2. They were first discussed 
in the papers [TT|.p4]. 

For Nf = 2Nc an exact type A (p^ — )■ — oo) solution of eq. (I3.25P is: 

P = N,+ ^N^ + Ql Q = Q^^m,^^^, 0<e<4. (3.26) 

Another solution with a qualitatively different UV behavior is: 

P = ^ + ce4^/3, c>0, Q = ±^. (3.27) 

One can check that in these solutions the radial coordinate moves all over the real axis 
and that for p — )■ — oo the solutions take the same form, but as anticipated above, differ 
substantially in the far UV, for p — )■ oo. Also, for the case Nf = 2Nc the papers |103] 
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discuss some extra solutions apart from the ones mentioned, including, interestingly, the 
generalization to near-extremal solutionis. 

3.4.2 Asymptotic expansions of generic solutions 

Other solutions of interest have been discussed in the papers [TT], [95], [21]. We will summa- 
rize the results, but suggest to the interested reader to go over those papers for details of all 
the metric functions. 

In the UV (for p — )■ oo), two possible asymptotics were found, that were called Class I 
and Class II in [24J. Table [1] summarizes the situation. 

In the IR (p — )■ 0), three types of solutions were found, called Type I, II and III (there 
exist other, qualitatively different solutions reported in [99] ). The function P{p) in these 
cases is, 

P = -Nfp + + ^4p>3 _ 2clNfP,p' + ^cl Qp„2 + Nj^ p' + O (/) , (3.29) 

for Type I. For the Type II asymptotics, we assume that this behavior occurs when the IR 
is located at pm > po. Without loss of generality we can choose pm = 0. With this choice 
we then necessarily have po < 0. Expanding Q in fl3.22p around p = we obtain 



Q = bo + bip + 0{p'), (3.30) 



where 



bo = - coth(2p„) {Qo + - " ^ * ^ 



2 2 



h = - si^h2^(2p ) + - (2iVc - Nf) coth(2po) . (3.31) 



The metric for the simplest non-extremal solution can be written in terms of a constant ^ and a function 
J" = 1 - {zh/z)'^ as: 

ds'^ = z\ - J'dt'^ + dxl + dxl + dxl + Nci^J'^^dz'^ + l(d0'^ +sin^ 9d(p'^) + 

\z^ 4 

+ ^^{de'^ + sn?M'f>'^) + ^{dilj + cosedLp + cosM(pf^ . (3.28) 

The solution also contains non-trivial RR F(3) and dilaton, see [TT] for details. Different features of this black 
hole solution have been analysed in |57j . |102j . An important remark is that the theory is in a Hagedorn 
phase and, indeed, the temperature coincides with the Hagedorn temperature of Little String Theory. For 
this reason, this solution is a bit problematic for studying the effect of quarks in a field theory plasma, unlike 
the finite temperature solution of section 12.51 
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I 


II 


< 2N, 


P r^Q \2Nc-Nf\p 

e^'' ~ I (27V, - Nf) p 
- N^ 

2N^[2N,-Nf)p 

a - {2Nc - Nf) e-^^P-P^^p 




>2N, 


Pr^-Qr^ \2N^ - Nf\p 

y ^ 1 (AT^ _ AT,) 

„4f _ , e4(''-''°)sinh2(2po) 

2(iVc-Ary)2(Ar^-2Ar^)p 
a e-'^^P-P°^p 


P ~ c+e^''/3 

e^f ~ c+e^p/^ 
Y ~ |c+e4^/3 
e4(<A-</'o) ^ ]^ 

a ~ 2e-2(''-''°) 


= 2iV, 


P-N^ + VN! + Ql-^jBm 

e4{<^-^o) _ e4{p-p„) gi^i^2 ^2pj (y^ 





Table 1: The two classes of leading UV behaviors. 
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Looking for IR solutions of fl3.25p we find that we must require that > 0. The correspond- 
ing asymptotic solution then takes the form 

P = Q + h,p'/'-^{hl + 12bo{h + Nf))p 

OOq 

+ ^ (5/1? + 6(56i + 2Nf)bo - ml coth(2p„)) p3/2 + oip"), (3.32) 

where hi is an arbitrary constant. Note that this expansion for P admits a smooth limit 
when po — )■ — oo and so it is valid for both solutions of type A (p^ — )• — oo). 

Finally, for Type III asymptotics we consider pi^ > po and we take Pjr = 0. In terms of 
the expansion (13.301) this requires that Bq = 0. We then find, 

P = /.ipV3 - ^p - ^ coth(2p„)p^/^ - ^ (506? - 18iVj) p^/3 + Oip% (3.33) 

where /ii 7^ is an arbitrary constant. 

To leading order the solutions for large p - UV solutions - are quoted in table [H It is 
the presence of subleading terms that allow the smooth numerical interpolation with three 
possible IR behaviors discussed. 

The physics of the dual field theory encoded in these solutions was discussed in detail in 
the papers [H] , [21] , [95] by computing various observables using the string solution of eq. 
(13. ip evaluated on We move now to discuss general features of the 

dual field theory. 

3.5 The dual field theory 

The proposal here is the following: without the addition of the flavor branes, the field theory 
is known to be a twisted version of six-dimensional Yang-Mills, or as we discussed above, a 
four dimensional QFT with an infinite number of massive fields. See section 13.2.21 To get an 
intuitive understanding of the modifications of the dynamics produced by the "quark" fields 
(that feature below), we will consider that all the infinite massive fields are chiral multiplets 
and then argue that the dynamics is ruled by a lagrangian of the form: 

L = Tr[-^F^^,-iXrD^X+L{<!>,,W)]^ j j d''e<!>y<!>k+ J ■ ■ . 

(3.34) 

When the flavor branes come into play, we are adding "quark superfields" that are realized 
as the open strings going from the non-compact flavor branes to the compact (or wrapped) 

-"^^Finding a numerical interpolation between the IR solutions and the solutions of Class I in the far UV is 
(numerically) delicate. One can see some plots in section V of the paper |100j . 
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color branes (as usual, the open strings that begin and end on a flavor brane decouple and 
do not contribute to the four dimensional dynamics) . More concretely, we add the quark and 
anti-quark superflelds [Q, Q) and propose that we have a lagrangian for the massive flelds 
interacting with the quark-antiquark superflelds Q, Q schematically of the formic 

L($fc, W,Q,Q) = J2j d^'e^ie^^k + f^k J (fOQ^kQ + /x^l^fcl^ + ■ ■ ■ , (3.35) 

k 

and canonical kinetic terms for {Q,Q). In this system, the SU{Nf)L x SU{Nf)ji symmetry 
is explicitly broken to the diagonal SU (Nf) £> by the presence of the coupling Q^^Q. In this 
respect, the theory is qualitatively different from M = 1 SQCD. 

One may be interested in the theory at low energies and hence integrate out the massive 
flelds (either massive vectors or massive chirals) and after some algebra end with a theory 
of the form (again schematically): 

L = j d^e(Q^e^Q + g^e-^g) + j (feW^W + k j (fe{QQf , (3.36) 

where we have a (naively irrelevant) deformation of A/" = 1 SQCD. 

We emphasize that this is an intuitive way of understanding the fleld theory dual to the 
flavored system described above. As we will summarize below there are various observables 
that can be computed that match the predicted (or expected) result. So, the precise dual 
QFT should be something similar to what we described above, or at least with the same 
qualitative physics. 



3.6 Checks and predictions 

This subsection summarizes results developed in the papers [11], [2l], [95]. There is a point 
that should be emphasized here. All the solutions to the BPS equations or the master equa- 
tion f l3.25p that have been found up to the time of writing this review, present a singularity 
in the IR (typically at p = 0). In spite of this being a "good singularity" according to 
some criteria developed in the literature [55]; the presence of the singularity makes the in- 
terpretation of IR observables a bit unclear. In other words, though one gets the "correct or 
expected" result one should perhaps handle those particular computations with care. 

Let us then concentrate on various quantities computed in the UV and then we will specify 
some that are mostly influenced by the IR of the geometry. 



3.6.1 Beta function and anomalies 

The gauge coupling and the theta angle of the dual QFT can be deflned as explained in 
various places, see for example section 4.1 of the paper [95j or section 5 of the paper [24]. 

^^For more details, see the appendix A of the paper [21] 
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One gets, after some algebra, that the gauge couphng is related to the functions of the 
background as: 

^ = e-^P . (3.37) 

Choosing a particular radius-energy relation, that was discussed in |104] . |105] . one can 
compute the variation of the coupling with respect to energy. Using the solutions where 
the dilaton asymptotes to a linear function {e^'^ ~ ^) and working to leading order in an 
expansion in inverse powers of the radial coordinate, we get 

/33.2 =h2N,-Nf) , (3.38) 

that coincides with the result predicted by the NSVZ result, once we assign anomalous 
dimensions to the quark superfields IQ = 1q = ~\- 

Similarly, one can define a geometrical quantity that can be associated with the quartic 
coupling. See section 4.2 of the paper [95]. The beta function can be computed using the 
anomalous dimensions discussed above and again get matching with the interesting fact that 
for Nf > 2Nc the quartic coupling is irrelevant, for Nf < 2Nc the coupling is relevant while 
for Nf = 2Nc the coupling is not running. See [106] for a nice explanation of this fact. 

One can also assign a value of the R-charge to the quark superfields to get the correct 
R-symmetry transformation properties of the quartic superpotential of eq. (13. 36 p . that is 
R[Q] = R[Q] = |- This predicts that the R-symmetry anomaly, the triangle with one R- 
current and two gauge currents is proportional to the quantity {2Nc — Nf) times the phase 
by which we are rotating the fermions. This is the precise result that the string background 
gives. Indeed, if we compute the 6-angle as explained in section 4.1 of the paper [HS] or in 
section 5 of the paper [21], we will get 

e = ^^^(2iV,-iV^) , (3.39) 

where we associated with the change in phase of the fermions in the quark multiplet and 
the gauge multiplet to get perfect matching. In the same vein, it is possible to attempt a 't 
Hooft matching of anomalies, that is of triangles involving three global currents. The reader 
will find it quite instructive to go over section 4.7 of the paper [95]. There, a detailed study 
of the matching of the correlator of three global currents -some of them corresponding to 
discrete symmetries, some of them being continuous symmetries- is presented. The treatment 
is performed in the case of the Type A backgrounds, that are characterized by the fact that 
the functions a = 6 = in eq. (13. ip . This translates to the fact that the R-symmetry is 
broken to Z2Nc-Nf without the further (spontaneous) breaking to Z2. 
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3.6.2 Seiberg duality 



It is known tliat Seiberg duality manifests beautifully in a QFT like the one of eq. f l3.36p . 
This is explained in section 1.10 of |106] . The backgrounds discussed here show this in a 
very nice way. Indeed, as discussed for example in the paper [21] , we can see that the master 
equation (and the whole system) is invariant under the change, 

P ^ P, Q ^ -Q, -a, ^ Nf - N^, Nf Nf , (3.40) 

while all other functions are invariant. Geometrically, this change amounts to swapping the 
two S"^ in the background, namely those parameterized by 9, ip and 9, (p in (13. ip . This should 
be interpreted as follows: suppose that we are presented with a background, representing 
the dynamics of a field theory with N^. colors and Nf flavors. This implies that we have a 
particular solution to the master equation (I3.25P for the function P. With this solution, and 
applying the changes of eq. (I3.40p . we can construct another solution, that will be related 
to the first one by a differomorphism and that will describe the physics of a field theory 
with Nf — Nc colors and Nf flavors. Various aspects of this interesting duality have been 
discussed in the papers [95], [21], [72] and probably elsewhere. 

The implementation of Seiberg duality in A/" = 1 subcritical string models was discussed in 
[53], |107] . One can think of the sphere exchange mentioned above as being the geometrical 
version of the mechanism described in these papers. Interestingly, similar methods were used 
to propose a non-supersymmetric Seiberg duality in [108] . 

3.7 IR Physics: Domain Walls and some comments on Wilson/'t 
Hooft Loops 

One observable that can be computed and that strongly depends on the p — )• region of the 
solution (the IR) is the tension of domain walls. Indeed, domain walls can be thought of 
as D5-branes that wrap a three cycle inside the internal six dimensional manifold and that 
extend on two of the Minkowski directions (and time, of course). We can compute the tension 
of a wall by considering a probe D5-brane that sits on the manifold Eg = [t,xi,X2,9,ip,ilj], 
at constant p = 0, constant 9, p. The Born-Infeld action for this probe can be computed 
and one reads that the effective tension is given by (see for example section 5.6 of the paper 

my- 

which, when evaluated at p = 0, gives a constant proportional to 2N(. — Nf. This is a good 
example for an observable, since even when a singularity is present (this typically reflects in 
some of the functions of the background being divergent), the combination above is finite. 
This is typical of "good singularities". There are other observables that can be computed 
using D-branes or fundamental strings, examples of these are Wilson or 't Hooft loops. A 
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similar conspiracy of functions that avoids an infinite result occurs here. Nevertheless, one 
should be quite careful with these quantities as noted above. Indeed, it was found in section 
V of the paper [100] that for the particular case of the backgrounds studied in this section, 
the Nambu-Goto action for the fundamental string might cease to be a good approximation 
as the string develops a cusp when approaching the singularity. 

In other words, we believe that the solutions presented in this section surely capture 
correctly many UV aspects of the field theory, together with some IR observables. Probably, 
we could think of the presence of the singularity in the same way as we think about the 
singularity in the Klebanov-Tseytlin backgrouncO, that captures some of the physics, but 
some is lost and the singularity must be resolved. There are different ways of attempting a 
resolution of the singularity; for example, considering massive quarks. This is under present 
study. 



3.7.1 Wilson loops and first order phase transitions 

One can consider the situation (of course, this is an idealized situation) in which all the 
flavors are massive and with a fixed sharp mass Mq (corresponding to a given value of the 
radial coordinate, that we call p = mg). The way to model this in a first approximation is to 
consider Nf{p) = NfQ^p — mo), where 0(p) is a Heaviside step function. Once again, there is 
some dynamics that is being lost in doing this, for example, the matching of the derivatives 
of the solutions is not smooth at the point p = uiq and the curvature of the background 
is not well defined at that point. Nevertheless, it is possible (within this approximately 
correct way of proceeding) to find a solution that for energies below the scale set by mg 
corresponds to the theory without flavors, say those discussed around eq. (13.91) and far in 
the UV corresponds to the flavored theory, as represented by solutions of the Class I in 
table [H One can then compute the Wilson loop following the well-known prescription [109j . 
This was done explicitly in |110j . The qualitative result is the following: for a range of 
ratios between the mass of the quarks Mq and the value of the gaugino condensate set by 
the function a(p) one observes that the relation between the quark-antiquark potential Vqq 
and their separation Lqq presents a first order transitioij^. In other words, a point where 
is discontinuous. The same kind of behavior was observed in systems where a more 
careful study is possible. Indeed, in the backreacted Klebanov-Witten (see section |2]) and 
Klebanov-Strassler (see section [5]) models, it was possible to find the precise form for the 
function p(p) - that shows that the "Heaviside approximation" described above is not a bad 
one. The same qualitative behavior for the first order phase transition was found in the 
papers [25], [m], [TT2] . 

^^In contrast to what happens for the solutions discussed here, the Klebanov-Tseythn background presents 
a had singularity and the IR physics computed with that solution is not trustable. 

^^Notice that we should talk of a "Quantum" phase transition, as the system is at zero temperature. 
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This kind of first order transitions for the Wilson loop string configurations are by no 
means particular of systems with dynamical fundamental fields. In fact, they were first 
found in a different context in [113] . where a nice connection between these Wilson loops 
computations and the Van der Waals gas (paradigm of the first order transition) was put 
forward (further discussions can be found for instance in [111] . |100] ). Different examples of 
such phase transitions in systems without flavors have been worked out in [114] . 

The "morale" seems to be the following: when we have a physical system that has two in- 
dependent scales (that is two scales that can be tuned independently, in the present example, 
the mass of the quarks Mq and the gaugino condensate A^) the first order phase transition 
for the quantity Vqq{Lqq) will be present. Of course, like in any other first order transition, 
it will happen that the discontinuity in the derivative will disappear for some ratio between 
the scales mentioned above. 



4 A dual to a (2+l)-dimensional A/" = 1 SQCD-like 
model 

In this section we will study gravity duals to minimal supersymmetric theories in 2+1 di- 
mensions. These backgrounds can be obtained by wrapping D5-branes along three-cycles 
of manifolds with G2 holonomy [1151 1116[ 1117] . The corresponding field theory dual is a 
(2-1-1 )-dimensional = 1 supersymmetric U{Nc) Yang-Mills theory with a level k Chern- 
Simons interaction. Such a theory coupled to an adjoint massive scalar field should arise 
on the domain walls separating the different vacua of pure = 1 super- Yang-Mills in 3+1 
dimensions. The corresponding unfiavored background was studied in ref. [116] . where it 
was argued to be dual to a U{Nc) gauge theory with Chern-Simons level k = Nc/2. In what 
follows we will review a generalization of these results, following closely ref. [118j . We will 
present the deformation of the background of [1151 1116[ [TT7] induced by a smeared distribu- 
tion of massless flavors. In order to formulate these generalized backgrounds, let cr* and a;* 
{i = 1, 2, 3) be two sets of SU(2) left-invariant one forms, obeying: 

da' = eijk (T^ A (x^ , du' = eijk A . (4.1) 

The forms a* and a;* parameterize two three-spheres. In the geometries we will be dealing 
with, these spheres are fibered by a one- form A\ The corresponding ten-dimensional metric 
of the type IIB theory in the Einstein frame is given by: 



(is = 62 



dxl, + dr'+^—{a'f+^—{u'-A'f 



(4.2) 
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where 0(r) is the dilaton of type IIB supergravity and g and h are functions of the radial 
variable r. In addition, the one-form will be taken as: 

A' = , (4.3) 

with w{r) being a new function of r. For convenience in this section we will take Qs = a' = 1, 
as we did in section 121 The backgrounds considered here are also endowed with a RR three- 
form F3. We will represent F3 as the sum of two contributions: 

= -^3 + /s , (4.4) 

where (iJ-3 = and /s is the part of the RR three-form which is responsible for the violation 
of the Bianchi identity {df^ 7^ 0) and which is sourced by the flavor D5-branes. Let us first 
parametrize the component J-3 as: 

^ = -hu^ - B^) A {uj^ - B^) A {uj^ - B^) + A {J - B') + H , (4.5) 

where B^ is a new one-form and are the components of its field strength, given by: 

= dB' + ^e,,,B^ AB'' . (4.6) 

In ( 14. 5p H is a. three-form that is determined by imposing the Bianchi identity for J-3, namely: 

dJ^3 = . (4.7) 

By using (14.11) one can easily check from the explicit expression written in (14.51) that, in 
order to fulfill (14. 7p . the three-form H must satisfy the equation: 

dH =\f' AF' . (4.8) 
4 ^ ' 

In what follows we shall adopt the following ansatz for B^: 

B' = il-llrl a' , (4.9) 

where 7(r) is a new function. After plugging the ansatz of B^ written in (14. 9 p into (14. 6p . one 
gets the expression for in terms of 7(r): 

F' = ^drAa' + e,,, A a'^ , (4.10) 

where the prime denotes the derivative with respect to the radial variable r. Using this result 
for in (14. 8 p one can easily determine the three-form H in terms of 7. Let us parameterize 
H as: 

H = ^ ^ n{r) ei,k A A a'' . (4.11) 
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Then, by solving (14. 8 p for H, one can verify that 'H(r) is the following function of the radial 
variable: 

7/ = 27=^-67 + , (4.12) 

with K being an integration constant. 

Let us now consider the contribution to the RR three-form F3. As explained above, this 
contribution violates the Bianchi identity and is non-zero when flavor branes are present. 
Indeed, let us write the WZ term of the action of a system of flavor D5-branes as: 

Sji^vor = n f QAC,, (4.13) 

JMio 

with Q being a four-form with components along the space transverse to the worldvolume 
of the branes. Then, the coupling to the RR potential Cq written in (I4.13P gives rise to the 
following modified Bianchi identity: 

dFs = dfs = An'^n . (4.14) 

To write an specific ansatz for Q and we have to select some family of supersymmetric 
embeddings for the flavor branes. As explained above, this can be done by using kappa 
symmetry. In the simplest case one looks for massless embeddings, which extend along the 
full range of the radial coordinate r. Those are the configurations considered in [118j . in 
which the D5-brane is extended along the three Minkowski directions x'^, as well as along a 
three-dimensional cylinder spanned by r and two other angular directions. Actually, it was 
shown in [118] that these two angular directions could be the ones corresponding to and 
u^. The corresponding transverse volume for this configuration is just: 

Vol( yl'^ ) = A A A . (4.15) 

However, there is nothing special in our background about these directions. Indeed, both 
in the metric and in the RR three-form J-3, we are adopting a round ansatz which does not 
distinguish among the directions of the two three-spheres. Thus we could as well consider 
supersymmetric cylinder embeddings that span the 1, 1 or 2, 2 directions. The volume forms 
of the spaces transverse to these embeddings are clearly: 

Vol( yl'^ ) = Aa^ Aoo^ Auj\ Vol( 3^4'^ ) = A A 00^ A . (4.16) 

To construct a backreacted supergravity solution with the same type of ansatz as in (14. 2 p we 
should consider a brane configuration that combines these three possible types of embeddings 
in an isotropic way. The corresponding transverse volume form Vol( 3^4) of this three-branch 
brane system would be just the sum of the three four-forms written in eqs. (I4.15P and (I4.16p . 
The corresponding smearing form Q is obtained by multiplying by the suitable normalization 
factor, namely : 

^ = -7^^01(3^4) = -^e,,keuma^ Aa'' aJ Aoo^ , (4.17) 



63 



where the minus sign has its origin in the different orientation (required in the kappa sym- 
metry analysis of |118] ) of the D5-brane worldvolume with respect to the ten-dimensional 
space. It is now straightforward to use the Q written in f l4.17p and get an expression of /s 
whose modified Bianchi identity is the one of fl4.14p . One has: 

h = ^e,,,{u' - ^) A A a' . (4.18) 

Eq. f l4.18p completes our ansatz for the general flavored case. Using these expressions of the 
metric and RR three-form in the supersymmetry variations of the dilatino and gravitino of 
type IIB supergravity, after imposing that the background preserves two supersymmetries, we 
arrive at a system of first-order BPS equations. These equations, which are rather involved, 
have been derived and analyzed in detail in |118j . They admit several consistent truncations 
which lead to simpler solutions. One can, for example, first consider the unfiavored case 
Nf = 0. If, in addition, we require that the function g is constant and that the fibering 
functions w and 7 are equal, our ansatz reduces to the one considered in |115l I116[ I117j . 
Actually, in this case the BPS equations fix the value of g to be e^^ = Nc and, in order 
to have a regular solution one should take the constant k of (14.121) to be equal to 1/2 (k 
is related to the Chern-Simons level k of the dual field theory). Other unfiavored solutions 
exist and have been studied in detail in ref. [118] . Here we will concentrate on reviewing 
the case in which Nf 7^ 0, starting from a particular truncation of the BPS system which is 
very interesting and serves to classify the different more involved solutions in the UV. 



4.1 The truncated system 

In this section we will analyze the truncation of the general system of BPS equations that 
corresponds to taking w = 7 = k = 0. In this case the BPS equations of |118j for the 
remaining functions h and g of the metric and for the dilaton consistently reduce to the 
following simple system of differential equations: 

4,' = JV^e-"" - 5(jv^_4iV^)e-»-2'> , 

g' = - ^e^'-^'^ - N,e-''^ + ~ e'^-^^ . (4.19) 

By inspecting the system fl4.19p one readily realizes that there is a special solution for which 
the metric functions h and g are constant. Actually this solution only exists when < 2Nf 
and the corresponding expressions for g and h are the following: 

1 {ANf-N,''^ 



e 



= ANf - AT, , e'' = - , {N,<2N/), (4.20) 
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while the dilaton grows hnearly with the holographic coordinate r, namely: 

, _ 2{3Nf - 



[4Nf-N,]-2 



■r r + 



(4.21) 



Let us next consider solutions for which the function h is not constant. In this case we 
can use p = e^^ as a radial variable and one can define a new function F{p) as F{p) = e^^ . 
It follows from f l4.19p that the BPS equation for F{p) is now: 



dF (^-A^c)(2-{)-^F 
dp F + N,- 4Nf 



(4.22) 



while the equation for the dilaton as a function of p can be written as: 

d(f) _ N, 
dp 



F{F + N^- AN^ 



f) 



3 / 4:Nf 

1 - — ( 1 - ^ 
4p 



(4.23) 



Moreover, from the second equation in fl4.19p we can obtain the relation between the two 
radial variables r and p, namely: 



dr 



dp F{p) + iVe - ANf 



(4.24) 



Notice that the sign of the right-hand side of (14.241) could be negative when Nj ^ 0. This 
means that we have to be careful in identifying the UV and IR domains in terms of the new 
radial variable p. We can use the result of integrating eqs. f l4.22p - fl4.24p to obtain the metric 
in terms of the new variable p, which takes the form: 



ds"' 
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f dr 



A' 



(4.25) 



Let us now study the different solutions of eqs. (14.220 - 04.230 . 



4.1.1 Linear dilaton backgrounds 

When Nf = 0, eq. f l4.22p can be simply solved by taking F = Nc- However, it is clear from 
f l4.22p that in the flavored case F = Nc is no longer a solution of the equations. Nevertheless, 
there are solutions for which this constant value of F is reached asymptotically when p — )■ oo. 
Indeed, one can check this fact by solving f l4.22p as an expansion in powers of 1/p. One gets: 

F = Nc + NcNf- - ^NcNf{Nc-4Nf)\ + ■■■ , (p ^ oo) . (4.26) 
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By plugging the expansion fl4.26p into (14.231) one can prove that, when 7^ 2A^j, these 
solutions have a dilaton that depends linearly on p in the UV and, actually, one can verify 
that: 

^ + ■ ■ • , (p ^ 00) . (4.27) 



dp 2{N^-2Nf) 8{N,-2Nfy p 

Notice the different large p behavior of the dilaton in the two cases Nc > 2Nf and Nf. < 2Nf. 
Indeed, when Nc > 2Nf the dilaton grows linearly with the holographic coordinate p (the 
behavior expected for a confining theory in the UV) , while for Nc < 2Nf the field decreases 
linearly with p. This seems to suggest that the sign of the beta function of the dual gauge 
theory depends on Nc and Nf through the combination Nc — 2Nf. Actually one can verify 
by means of a probe calculation in the complete system that the beta function is positive 
for Nc > 2Nf and changes its sign when Nc < 2Nf |118] . 

Eq. (14.221) can be solved numerically by imposing the behavior (I4.26P for large p. Once 
F(p) is known one can obtain the dilaton 0(p) by direct integration of the right-hand side 
of ( 14.23p . The result of this numerical calculation was analyzed in detail in |118j . Let us 
only mention here that, in the most interesting case Nc > ^Nf, the function F diverges for 
p — )■ 0, while the dilaton (p remains finite for small p. This bad IR behavior of F is cured in 
the untruncated solution with the same leading UV form of F and (p but with ly, 7 7^ (see 
below). 

4.1.2 Flavored G2 cone 

Let us now consider the solution of the equations (14.221) and (14.231) that leads to a metric 
which is asymptotically a G'2-cone with constant dilaton in the UV. It can be checked that 
there exists a solution of (I4.22p which can be expanded for large values of p as: 

4 , , 15N^ - 39NcNf + 24N^f 

F = - p + A{Nf - Nc) + ^ + ■ ■ • . (4.28) 

3 p 



The corresponding expansion for 0(p) is: 

^iV,(]V, + 2iV^)4 + ... , (4,29) 



4 p 32 ^ ^ p2 



where is the constant limiting value of (p in the UV. In order to explore the asymptotic 
form of the metric for large p, it is convenient to perform a change in the radial variable, 
namely: 

P = It\ (4.30) 
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in terms of which the metric asymptotically becomes the one corresponding to the direct 
product of (2+l)-dimensional Minkowski space and a seven- dimensional cone with G2 holon- 
omy, namely: 



dxl, + {drf + ^ + Y (^^ - y ) • (4-31) 



To find the solution in the whole range of the radial coordinate one can numerically integrate 
the system fl4.22p - fl4.23p by imposing the asymptotic behavior f l4.28p to the function F{p). 
For Nc > 2Nf one can show that F{p) is well-defined for p > 0, while it diverges for p — 
(see |118j for further details). Notice that, at least in the unfiavored case Nf = 0, it is 
natural to regard these solutions with finite dilaton in the UV as corresponding to D5- 
branes wrapped on a three-cycle of a G2 cone, in which the near horizon limit has not been 
taken and, thus, as we move towards the large p region the effect of the branes on the metric 
becomes asymptotically negligible and we recover the geometry of the G2 cone where the 
branes are wrapped. 



4.2 The complete system 

Let us now consider the solutions of the BPS equations for our general ansatz. These 
complete BPS equations have been derived in the appendix A of |118] . Here we will restrict 
ourselves from now on to the cases with A^^^ > 2A^j, which are the ones that lead to more 
sensible solutions. As in the truncated case of section 14.11 we will use p = e'^^ as radial 
variable and F = e^^ as a function of p. In order to solve the general BPS equations we must 
impose initial conditions to the functions w{p) and 7(p) introduced in (14. 3 p and (14. 9p . and 
we must fix the value of the constant k of (14.121) . These initial conditions are determined by 
imposing some regularity requirements at p = that we now review (see [118j for additional 
details). First of all, we will demand that the function F approaches a constant finite value 
when p — )■ (z.e. F ~ Fo for p — > 0). In order to fix the value of the function w{p) at p = 
let us recall (see (14.30 ) that w parameterizes the one-form which, in turn, determines the 
mixing of the two three-spheres in the ten-dimensional fibered geometry. The curvature of 
the gauge connection (defined as in (14.61) with — )■ A^) determines the non-triviality 
of this mixing. When this curvature vanishes one can choose a new set of three one-forms 
in which the two three-spheres are disentangled in a manifest way and one can factorize 
the directions parallel and orthogonal to the color brane worldvolume in a well-defined way. 
From the wrapped brane origin of our solutions, one naturally expects such an un-mixing of 
the two S'^^s to occur in the IR limit p = of the metric. Moreover, by a direct calculation 
using (14.11) it is easy to verify that for w = 1 the curvature of the one-form A^ vanishes . 
Thus, it follows that the natural initial condition for w{p) is: 

w(p = 0) = l. (4.32) 
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Actually, the three-cycle that the color branes wrap can be identified with the one that 
shrinks when p — 0, which is the one given by: 

j: = {uj' = a'}. (4.33) 

In order to have a non-singular flux at the origin, the RR three-form F3 should vanish on S 
when p — >■ 0. It is easy to check that this occurs if the constant k takes the value: 



2 2N, 



(AM) 



c 



Actually, f l4.34p is also a necessary condition to have a flnite dilaton at p = 0. Indeed, it was 
shown in |118j that, in addition to (14. 34 p . the dilaton remains flnite in the IR if the function 
7(p) takes the following value for p = 0: 

7(p = 0) = 1 - ^ . (4.35) 

Eqs. (I4.32P and (14.350 provide the initial conditions for the functions w and 7 we were 
looking for. 

4.2.1 Asymptotic linear dilaton 

As explained above, we are interested in solutions of the BPS equations such that asymp- 
totically F is constant. Actually, by solving the BPS system in powers of 1/p, one can check 
that there are solutions in which F has the following asymptotic behavior: 

F = N, + ^ + ^ + ^ + --- , (4.36) 

p p^ p-^ 

where the coefficients oi, 02 and 03 are given by: 

ai = NfN, , 

a2 = -^N,Nf{N, -4Nf) , 



NfN, r 



21N^ - USNfNc + 240 A^? 



(4.37) 



Notice that the first two terms in (I4.36P and (I4.37P coincide with the one written in (I4.26P 
for the truncated system. Similarly, the functions w and 7 can be represented as: 



w 





b2 


+ 


h 


p 


p2 


p3 


^ + 


C2 


+ 


C3 








p 


P' 




p3 



7 = - + ^ + ^ + --- , (4.38) 
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where the coefficients bi and q are the following: 



fei = ci = - (iV, - 3Nf) , 



b2 = C2 = - (iV, - 3Nf) (iV, - 2Nf) 



A9N^ - ISAN^Nf + 20A Nf 



C3 



32 



{Nc -SNf) A9N^ - 208 NfNc + 252 



(4.39) 



Moreover, for Nc > 2Nf the dilaton grows linearly with p as in (14.271) . i.e. (f) ^ p/[2{Nc—2Nf)] 
for large p. 

The solution for the full range of the holographic coordinate can be found by numerical 
integration of the BPS system with the IR regularity conditions (I4.32p . (14.341) and (14.351) 
and with F{p = 0) = Fq finite. One has to perform an interpolation between the p — )■ and 
p — i- 00 behaviors by means of a shooting technique in which the only free parameter Fq is 
varied until a solution with F{p) ^ A^^^ for large p is obtained (which only occurs when Fq is 
fine tuned to a very precise value). 

After obtaining this solution of the equations of motion of the gravity plus brane system, we 
can see if it incorporates some of the features that the supergravity dual of 2+1 dimensional 
gauge theory plus flavors should exhibit. In particular, we can study the evolution of the 
gauge coupling constant with the holographic coordinate. In order to do that, let us consider 
a D5-brane probe extended along the three Minkowski directions and wrapping the internal 
three-cycle S defined in (I4.33P at a fixed value of the holographic coordinate p. By looking 
at the J-"^ terms in the DBI action of this probe, we get the value of the Yang-Mills coupling 
constant of the dual (2+l)-dimensional gauge theory, namely: 



9ym 



det(G3) d^^ 



P + 



:i-w) 



(4.40) 



where G3 is the induced metric on the three-cycle S and we have neglected all constant 
numerical factors. Due to our initial condition (14.321) . the right-hand side of (I4.40p vanishes 
for p = 0, which corresponds to having — 00 in the IR, as expected in a confining 
theory. Moreover, I/qym grows as we move towards the UV region p — )■ 00, in agreement 
with the expected property of asymptotic freedom. Other gauge theory observables for these 
backgrounds, such as the Wilson loops, can be also analyzed (see ref. |118j ). Notice that, 
despite the regularity conditions we have imposed, in the flavored case Nj ^ the explicit 
calculation of the scalar curvature for the linear dilaton solutions shows that the metric is 
singular at the origin of the radial coordinate. Notice that, as argued for other backgrounds. 
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it is physically reasonable to expect that massless flavors drastically alter the backreacted 
geometry in the deep IR. However, as our initial conditions are such that the dilaton is 
finite at the origin, the value of the gu component of the metric is also bounded and then, 
according to the criterium of [66], the singularity is "good" and the background can be used 
to extract non-perturbative information of the dual gauge theory. 



4.2.2 Asymptotic G2 cones 

When F{p = 0) takes values in a certain range, the solutions of the BPS equations lead to 
the metric ( I4.3ip at the UV, which is the direct product of 2+1 dimensional Minkowski space 
and a G2 cone. The solutions in this case are very similar in the UV to the ones discussed 
in subsection 14.1.21 (with better IR behavior) and we will not discuss them further here. Let 
us only mention that the asymptotic values of F, w and 7 for p — j- 00 can be determined 
analytically and are given by: 

F ^ ^p + A{Nj-N,) + ■■■ , 

w ~ ^ — — a + . . . , p ^ 00 , 

2p 

1 Nf 

+ (4,41) 



5 Flavors in the Klebanov-Strassler Model 

The so-called Klebanov-Strassler (KS) solution [12] is dual to a cascading, confining theory, 
and has been a popular and successful laboratory in which to study numerous issues related 
to gauge-gravity duahty and to cosmology. The gauge theory lives on a stack of regular and 
fractional D3-branes at the tip of the deformed conifold, as we now briefly review. 

The deformed conifold is a regular, six dimensional, non-compact manifold defined by the 
equation 2:1^2 — z^z^ = (x^ in C^. When the complex deformation parameter fi is turned 
off, it reduces to the singular conifold, which is invariant under complex rescaling of the Zi. 
The base of the conifold has SU{2) x SU{2) x f/(l) isometry and 5^ x topology. The 
deformation parameter breaks the scale invariance, produces a blown-up at the apex of 
the conifold and breaks the U{1) isometry to Z2. 

The low energy dynamics of regular and M fractional D3-branes on the deformed 
conifold is described by a cascading A/" = 1 4d gauge theory with gauge group SU {N + M) x 
SU{N) and bifundamental matter fields A, B transforming as SU{2) x SU{2) doublets and 
interacting with a quartic superpotential Wkw = AiBkAjBi. The dual to this theory 
is the KS solution [12], that is relevant for the N = nM case, where n is an integer. The 
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related theory develops a Seiberg duality cascade which stops after n — 1 steps when the 
gauge group is reduced to SU{2M) x SU{M). The regular KS solution precisely accounts 
for the physics of an A ^ i?-symmetric point in the baryonic branch of the latter theory, 
which exhibits confinement and — )■ 1'2N ^2 where the second breaking is due to 
the formation of a gluino condensate (AA) ~ A^^. The complex parameter e is the geometric 
counterpart of this condensate. 

In this section, we will discuss how the solution is modified when a smeared distribution of 
D7-branes is introduced. In the dual theory, they correspond to fundamental fields, but the 
precise way in which they couple to the rest of fields depends on the D7-brane embeddings, 
as we will discuss below. In what follows, we only discuss cases in which the flavor D7-branes 
do not break any supersymmetry, such that the four-dimensional = 1 of the KS solution 
is preserved. The material we summarize in this section was developed in [G^, [119] . 

5.1 Backreaction with non-chiral flavors 

5.1.1 Brane embeddings 

Let us start by choosing an appropriate family of supersymmetric D7-brane embeddings. 
A particularly interesting example is given by D7-branes wrapping the holomorphic 4-cycle 
defined by an equation of the form [73] : 

Zi-Z2=fX, (5.1) 

where /z is a constant. It was shown in [73] that this embedding is K-symmetric and hence 
preserves the four supercharges of the deformed conifold theory. 

A D7-brane wrapping the 4-cycle defined above is conjectured to add a massless (if ^ = 0) 
or massive (anti) fundamental flavor to a node of the KS model. The resulting gauge theory 
is said to be "non-chiral" because the flavor mass terms do not break the classical flavor 
symmetry of the massless theory. The related perturbative superpotential is, just as in the 
singular conifold case [H], which we wrote in fl2.14p . The complex mass parameter m in W 
is mapped to the geometrical parameter fi. The different flelds are summarized in the quiver 
diagram of flgure [61 

By acting on this flducial embedding (15. ip with the generators of the broken symmetries, 
we can build the family of embeddings over which we want to smear. This is the obvious 
generalization to the deformed conifold case of the discussion in section [2] and, in fact, a 
generic non-chiral embedding is still given by (I2.17p . 

5.1.2 The ansatz 

We now write the ansatz for the metric and forms. It is similar to the ansatz for the KS 
solution, but, due to the presence of D7-branes, the RR one-form is non-trivial and the 
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Figure 6: The quiver diagram of the gauge theory. Circles are gauge groups, squares are 
flavor groups, and arrows are bifundamental chiral superfields. Nfi and Nf2 sum up to Nf. 



dilaton runs. It is useful to introduce the gi one-forms used in [I2] : 



— sin 9i difi — cos ijj sin 62 d(p2 + sin ijj d62 2 '^^1 ~ ^ ^2 d(p2 — cos ip d02 

sin 9i difi + cos ip sin 6'2 dip2 — sin (i^2 4 c^^i + sin sin 62 dip2 + cos ip d92 
, 9 = 



9 



^/2 



g = dip + cos 61 difi + cos 62 dip2 
The Einstein frame metric ansatz i?^^'- 



(5.2) 



di' 



3/C3(r) 



{dT' + i9r) + 



+ cosh^ (0 ((/)^ + (gr) + sinh^ (0 {{g^ + (/)^) 



(5.3) 



where fi is the complex deformation parameter of the conifold, dxf^ denotes the four- 
dimensional Minkowski metric and }C(t), h{T) and the dilaton are unknown functions 
of the radial variable to be determinedr^ 



^^A more generic form of the ansatz was used in |64| . [111] . By requiring supersymmetry and performing 
some algebra, one ends up with (|5.3p . We wih skip those intermediate steps here for the sake of briefness. 

The relation of the Zi complex variables as used above to the T,9i,(f 1,62,^2, coordinates can be 
found, for instance in . The embedding equation (I2.17|) expressed in terms of the "deformed conifold r 
variable" looks the same in terms of the "backreacted ansatz r variable" . See [111] for details. 
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For the forms we will adopt the following ansatz: 



B2 

Fi 
F3 



dh ^(r) A dx^ A ■ ■ ■ A dx^ - h—e ^3 sinh^ r JC^gi A g2 A A A . 

16 



9sa- 



2 

9sa- 



fg'Ag^ + kg^Ag' 



drAifg'Ag' + k g' A g') + -{k - f) g' A {g' A g' + g' A g') 



Arc 



+ '-^^f)g^ A / + (1 - F + ^^^,),3 ^ ^4 



47r 



+ 



+FdTA{g'Ag' + g^Ag')}, 



(5.4) 



where / = /(t), A; = A:(t), F = F{t) are functions of the radial coordinate (and where the 
dot denotes derivatives with respect to r). We have implemented the self-duality condition 
for F5. 

Notice that, consistently, dFi = —gs^, where Q is the symmetry preserving D7-brane 
density distribution form analogous to (I2.3ip : 



Q = ^ ^p(r) (sin 9id9i A dipi + sin 92d92 A dip2) - p(r)rfr A g^ . 



(5.5) 



When quarks are massless |64j, one just has p(r) = 1, whereas p(r) becomes non-trivial 
when quarks are massive. In what follows, we will keep p(r) generic. We refer the reader 
to for the computation of p(r) from the massive non-chiral brane embeddings (I2.17p . 
The source contributions to the modified Bianchi identities for F3 and F5 



dF^ = ifg A Fi - g^n A B2 , 

dFs = H3AF3 -\gsnAB2A B2 



(5.6) 



follow from the WZ term of the smeared D7-brane action [111] . Given (15. 4p and (15. 5p . the 
equations (15. 6p are satisfied provided: 



h^e-'-T sinhV/C^ = const - -{gsa'Mf 
16 4^ ^ 



f-{f-k)F + ^p{r)fk 
Air 



(5.7) 



5.1.3 The BPS equations 

By requiring the vanishing of the bulk fermionic supersymmetry variations, one finds a set 
of first order BPS equations. The computation is lengthy but straightforward and was 
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carried out in [6l] (since in that paper p(r) = 1, the substitution Nf — )■ Nfp{T) has to be 
implemented in the equations of [M])- In the present notation, the differential equations are: 

4:71 

k = e^^F+^^^^/Voth^-, 
/ = e^(l-F+^£^fc)tanh^I, 
F = le~%k-f), 

§ = +^-cothr , (5.8) 

/C sinh r 3 

supplemented by the algebraic constraint 

e-^(k -f)= tanh - - 2F coth r + ^'^^^^^^ \k tanh - - / coth -1 . (5.9) 

Quite remarkably, the equations fl5.8l) - fl5.9p can be (almost) explicitly integrated. In the 
following, we will use notations similar to those employed in section [2j We introduce an 
arbitrary value of the radial coordinate t^, at which the dilaton is 0=,,. Then, we can write 
the dilaton as: 

where we have introduced the deformation parameter which weighs the flavor loops as: 



Let us also introduce a function: 



A* with K = A7rgsMe'^' . (5.11) 



r^(r) = e.e'^-'^* [ (sinh2e - 20^(0^^^ • (5.12) 
Then, we can integrate for the rest of the functions of the ansatz: 

^ [ sinh 2r - 2r + r7(r)] ^ ^ sinhr-r 

2^ sinhr ' 2 sinh r 

J T coth r — 1 , , , ,T coth r — 1 , 

/ = coshr-1 , k = coshr + 1 . 5.13 

2 smh T 2 smh r 

Finally, the function h can be obtained by integrating (15. 7p . The KS solution without flavors 
is obtained by taking = 0, such that the dilaton is constant and rjlr) = 0. For p(r) = 1, 
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we find the solution backreacted with massless flavors [6^ . In this case, the integrals for the 
dilaton and //(r) can be explicitly performed: 

ri(T) = e,e'^-**(sinhV-r2) , 

e*-^* = -— ^ r. (forp(r) = l) (5.14) 

1 + e*(r* — T) 

In this massless case, the solution has a curvature singularity in the IR r = 0. Some cases 
where p(r) is non-trivial were discussed in 

5.1.4 Some physical features 

The solution presented in the preceding sections has been used to extract some of the physics 
encoded in the unquenched background. In [6lj, the running of the couplings and anomalies 
were discussed. As anticipated above, in the solution with massive flavors was found. 

Quark masses erase the IR singularity in the same way as explained in section 11.51 or in 
section 12.3.21 Quark-antiquark potentials, screening lengths and associated quantum phase 
transitions were discussed in the same paper. Finally, in [70], it was computed how the 
screening effects due to unquenched fundamental matter affect the mass spectra of the KS 
model, with results similar to section 12.41 Due to space constraints, we cannot go explicitly 
through all of these features and we refer the interested reader to the original papers. Here, 
we will just briefly discuss how the solution captures the phenomenon of a duality wall [M] 
and how gauge groups ranks change upon Seiberg duality. 

We will make use of the following holographic formulae, which can be derived in the 
M = 2 orbifold case by looking at the lagrangian of the low energy fleld theory living on 
probe (fractional) D3-branes: 



An^a' Js2 ^ 2 



gYhSY ' gvuisY gs ' gvivSY gvuisY gs lAn^a' Js2 ^ 2 J ' ^^'^^^ 

The labels (l), (s) in f l5.15p refer to the gauge group with the larger or smaller rank. Strictly 
speaking, these formulae need to be corrected for small values of the gauge couplings and 
are only valid in the large 't Hooft coupling regime (see |120| \10Q\ l63]). which is the case 
under consideration. Moreover, they are also expected to be precise just in the UV region, 
where the cascade takes place and the region on which we will focus below. The expressions 
f l5.15p give positive squared couplings only if the expression inside the square bracket is in 
the range [— |, |]. Deflne: 

Mr) [ B, = ^/ = ^ , 6o = &o - N e [0, 1] , (5.16) 

where [6o] denotes the integer part of h^. In order to get the explicit expression for fog we have 
integrated over the S'^ parameterized by 6i = 02, ipi = 27i — ip2, ip = const [6l], considered 
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the UV limit of flS.lOp . fl5.13p such that p(r) ^ 1 and f k e'^ij — l)/2 and inserted the 
definitions (15.1 ip . Now we see that what we have to insert in (15.150 is indeed 6o- This is 
the physical content of the cascade: at a given energy scale we must perform a large gauge 
transformation on B2 in supergravity to shift / i?2 by a multiple of An'^a' to get a field theory 
description with positive squared couplings. 

Let us restrict our attention to an energy range, between two subsequent Seiberg dualities, 
where a field theory description in terms of specific ranks holds. When fiowing towards the 
IR, bo decreases from 1 to 0. From (I5.15p and inserting the solution, we can find an expression 
for each of the gauge couplings: 

V = ^(1 - - r*)) 60 , ^ = ^(1 - - r,)) (1 - bo) . (5.17) 
A/ A* As A=K 

In this energy range, the coupling A/ starts different from zero and fiows to 00 at the end 
of this range, where a Seiberg duality on its gauge group is needed. The coupling of the 
gauge group with smaller rank is the one which starts very large (actually divergent) after 
the previous Seiberg duality on its gauge group, and then flows toward weak coupling. 

The qualitative picture of the RG flow in the UV can be extracted from our supergravity 
solution even without discussing the precise radius-energy relation, simply recalling that the 
radius must be a monotonic function of the energy scale. First, notice that at a flnite r (and 
therefore at a flnite energy scale Ejjy)^ the dilaton diverges making both gauge couplings 
diverge. This happens at: 

Tdw = r^^ • (5.18) 

e* 

From fl5.16p . we see that the derivative ^ grows unbounded near t^w^ meaning that the 
interval (in r) between Seiberg dualities becomes shorter and shorter. The Seiberg dualities 
pile up the more we approach the UV cut-off Euv The picture which stems from the flavored 
Klebanov-Tseytlin/Strassler solution is that Tdw is a so-called "Duality Wall", namely an 
accumulation point of energy scales at which a Seiberg duality is required in order to have a 
weakly coupled description of the gauge theory |121j . Above the duality wall, Seiberg duality 
does not proceed and a weakly coupled dual description of the fleld theory is not known. 
See Figure [3 Nevertheless, in full analogy with the discussion of section |2l the derivative of 
the holographic a-function changes sign at a flnite distance in r below Tdw so one should not 
trust the solution all the way up to the singular point Tdw 

Duality walls were studied in the context of quiver gauge theories flrst by Fiol |122j and 
later in a series of papers by Hanany and collaborators [123]. To our knowledge, the solution 
above is the only explicit realization of this exotic ultraviolet phenomenon on the supergravity 
side of the gauge/gravity correspondence. 

To end this section, we discuss how the effective number of regular and fractional D3- 
branes change when undergoing a step of the cascade of Seiberg dualities. We will not 
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Figure 7: Qualitative plot of the running gauge couplings as functions of the logarithm of 
the energy scale in the cascading gauge theory. The blue lines are the inverse squared gauge 
couplings, while the red line is their sum. 



compute the explicit shift in r but rather the shift in the function / k). From fl5.16p . we 
have: 



&o(r)^&o(r') = &o(r)-l 



fir) ^ fir') = fir) 



TT 



gsM 



(5.19) 



On the other hand, we compute the effective number of branes at a given energy scale by 
integrating the appropriate RR-forms: 



1 



^^eff{r) 
Me f fir) 



^ ' F3 = M 



iVo + 



gsM^ 



TC 



f + 



9sNf 2 
An ■' 



A-n'^Qsa' J S3 



1 + ^/ 

271 . 



(5.20) 



In these expressions we have substituted (15.41) . (15.71) and already taken the UV limit f = k 
and p{r) = 1. The for the second integral is the one parameterized by 62 = constant, 
(p2 = constant. Notice that N^jf and M^jf are not quantized. This is because they are 
Maxwell charges, as opposed to Page charges. See [M] for thorough explanations. 

We can compute how N^ff, M^jf vary in a Seiberg duality step (I5.19p . A bit of algebra 
shows that: 



M. 



eff\ 



M, 



T 



M, 



eff 



r 



2 ' 



N^ffir) ^ N,^^(r')=N,fjir)-M,jf{r) + ^ , 

whereas Nj remains unchanged. A careful analysis in [61] showed that this is in full agree- 
ment with field theory expectations. 
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5.2 Backreaction with chiral flavors 



In a remarkable paper [119] . Benini discussed the solution dual to having smeared chiral 
flavors on the conifold. In the probe approximation, the D7-brane embeddings that cor- 
respond to chiral flavors were discussed in [H], |124] . How these flavors transform under 
the gauge groups is shown in the quiver diagram I5.2[ In this case, the quiver theory is not 
self-similar under the duality cascade; in each step of the cascade a meson fleld is generated. 
Its couplings to the rest of the flelds are, however, irrelevant |119] . 




Figure 8: Quiver diagram of the KS theory with chiral flavors. 

The backreacted solution of jllQj uses the singular conifold and therefore it can be con- 
sidered as the deformation of the Klebanov-Tseytlin solution [125] due to smeared chiral 
flavors. From the gravity point of view, the extra complication with respect to section 15.11 
is that the worldvolume gauge fleld on the D7s has to be turned on. In fact, this is crucial 
when matching the shifts in the ranks of the gauge groups upon Seiberg dualities to the 
supergravity background (there are subtle differences with respect to the non-chiral case). 
We will not report further on this solution here, but refer the reader to |119] . 

6 Models with cohomogeneity 2 

In this section we present some situations in which, even smearing the flavor branes, the 
system cannot be reduced to a one-dimensional problem. In fact, the different flelds will de- 
pend on two different radial coordinates and, accordingly, one has to solve partial differential 
equations rather than ordinary differential equations. 

In order to provide a heuristic picture, the situation is depicted in flgureO Concretely, we 
will refer here to the model of section 16. H but the situation is very similar for all the cases 
discussed in this section. 

In flgure [9l the color branes are placed at the tip of a Calabi-Yau (a is a radial coordinate 
along the CY, the rest of the directions of the CY are omitted from the plot). The p — 02 
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Figure 9: A qualitative plot of the situation with cohomogeneity 2 models. The red dot 
in the center represents the color branes and each vertical line is a flavor brane. Taking a 
lot of them smeared along the 02 angle, the rotational symmetry associated to this angle is 
effectively recovered. Any function of the ansatz depends on the radial coordinates p and cr. 

plane is transverse both to the color branes and to the CY. Each flavor brane lies at a 
point in this plane and is extended along a. Distributing the flavor branes along 02, it is 
possible to recover (in the smeared limit) the associated U{1) isometry. On the contrary, as 
is apparent from figure [H there is no way in which one can place the flavor brane to recover 
the full radial symmetry. Hence, the solution associated to this brane configuration must be 
cohomogeneity two, meaning that all functions of the eventual ansatz will depend on p and 
a. 

As a matter of fact, if one wishes to construct a deformation of AdS^ x 5*^ with smeared 
flavor such that the supersymmetry preserved is A/" = 2 (rather than A/" = 1 as in section |2]) , 
the solution would have cohomogeneity two and, presumably, would share some similarities 
with the examples presented in this section. This is an interesting open problem for the 
future. 

6.1 A dual to (3+l)-dimensional TV = 2 SQCD-like theory 

In this section, we study the dual solution to the brane intersection summarized in table El 
The gauge theory lives on N^. D5-branes wrapping a two-sphere with the appropriate twisting 
to preserve eight supercharges, i.e. A/" = 2 in the effective four- dimensional low energy theory. 
Geometrically, it corresponds to wrapping the branes along a compact SLag two-cycle inside 
a non-compact Calabi-Yau two-fold. This leaves two flat transverse dimensions which are 
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identified with the moduh space corresponding to giving vevs to the complex scalar inside the 
M = 2 vector multiplet. The Nf flavor D5-branes do not further break supersymmetry and 
provide fundamental hypermultiplets in order to build M = 2 SQCD. They are extended 
in the non-compact a direction and, thus, their volume is inflnite, making exactly zero 
the effective four-dimensional gauge coupling living on them. They would provide a global 
symmetry group U{Nf) if they were placed on top of each other, but due to the smearing, 
only U{l)^f is left. The dual solution without flavors was found in |126] . and the flavored 
case was discussed in |127] . 



CY2 M2 





Xl,3 


a 


01 


9 




P 


02 


D5 


















Nf D5 


















Table 2: A scheme of the setup: for the brane conflguration, a line — means that the brane 
spans a non-compact dimension, a point • that it is point-like in that direction, a circle Q 
that it wraps a compact cycle and ~ indicates smearing in the direction. Above, it is shown 
which directions spanned the Calabi-Yau and which the transverse plane before backreaction. 



We start by writing an ansatz for the metric consistent with the symmetries of the problem. 
In Einstein frame: 



ds^n = gsN^a'e^ 



^ 7^^1,3 + ^(de'^ + sin^ dd(p^) + 



2$ 

(6.1) 



— 2$ 

+e-'*(rfp2 + p'dcg) + [da^ + a\d<p^ + cos^^)^ 



where z and $ depend on both radial coordinates p, a. The Calabi-Yau twofold directions 
are 0<^^<7r, 0<(^< 27r, < a < 00, < 01 < 27r (of course, in this solution with fluxes 
there is not a Calabi-Yau any more, but it can be thought of as a deformation of the Calabi- 
Yau that was present before backreaction). The coordinates < p < 00, < 02 < 27r span 
the transverse two-dimensional plane, so they should be identifled with the moduh space, 
and therefore rotations in 02 are related to the U{1)fi symmetry of the fleld theory. Out of 
the SU{2)r symmetry, only its diagonal U{l)j is manifest in the geometry, as rotations in 
01. The extra 50(3) isometry which acts on 9, (p, 0i does not play a role in the low energy 
N = 2 SQCD theory [T26] . 

As anticipated in tabled we want to consider a set of Nf D5-branes extended in xi^s, a 
and wrapped in 0i. They lie at flxed p = pg, where pq is proportional to the modulus of 
the mass of the fundamental hypermultiplets. These D5-branes are homogeneously smeared 
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over the S"^ parameterized by 6, ip and on the angle 02, which corresponds to the phase of 
the mass of the hypers. This distribution is described by the four-form: 



n = ^5{p - Pq) sin edp A rf02 A dO A dip , (6.2) 



6TT 



such that the source-modified Bianchi identity for F(3) reads: 

(iF(3) = 2k(io)T5(^ = g,a' (p - Pq) sin 9dp A A dO A d(p . (6.3) 
We can write an ansatz for F(3) consistent with this expression: 

F(3) = NcQsOi' —g'd(j)2 A dp A {dcpi + cosOdip) — gd(f)2 A da A {dcpi + cos9d(p) + 

(6.4) 



Nf ^ ^ 

+ (9 + 771^0(P - Pq)) sin ^#2 AdOAd^ 



2N, 

where B is the Heaviside step functioij^. g a new function of p and a that needs to be 
determined, and we have introduced the following notation for the partial derivatives: 

' = dp, - = 8^ (6.5) 

The next step is to insert the ansatz (16. ip . (16.40 into the type IIB supersymmetry transfor- 
mations Sipfi = SX = 0, as outlined in section 11.4.11 This procedure was carefully performed 
in [127] . whereas here we just quote the resulting system of first order equations: 

Nf 

9 + ^0(P - Pq) = -pz' , g' = -2e-^^pa^ , 

e^* = ^ , g = -z-'e-'%{g + ^Q{p - pq)) + 2z'' pae-'"^' . (6.6) 

It is easy to check that the last equation is not independent of the previous ones and that 
equations (16. 6p ensure the equation of motion for the 3-form (i(e* *-F(3)) = 0. This system 
of equations can be recast as a single, non-linear, second order PDE for z(p, cr): 

Nf 

^2N~^^^ - Pq) + P^i^ - ^'^) = ^(P^^ + + pz") ■ (6-7) 

Once z{p, a) is computed, g and $ are read from (16.60 . In general, the equation (16. 7p cannot 
be solved explicitly. In the unflavored case Nf = 0, there is in fact an exact solution [126j 
(see |127j and the first paper of |105] for the adaptation of the solution [126j to the present 
coordinate system). 



^"'^Notice that, as opposed to section [3.7.11 where a Heaviside function was introduced as an approximation 
to the effect of the massive flavors, the Q here is exactly what comes from the family of D-brane embeddings 
considered, since they all lie at fixed p = pq. 
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Equation (16. 7p . however, can be studied numerically |127j . We will not pursue that here, 
but we will verify using (16. 6p that the expected beta-function for the gauge coupling stems 
from the differential equations. In order to read the effective four-dimensional gauge coupling 
from the geometry, we consider a "color" D5-brane probing the Coulomb branch of the theory, 
namely, a D5 wrapping the S"^ parameterized by 9, (p, sitting at cr = [105] . After integrating 
the volume of the S"^, we find: 



9ym 



Nc 
47r2 



(6.8) 



Thus, in order to understand the running of the coupling it is not necessary to know the 
geometry everywhere, but just at cr = 0. From the second equation of (16. 6p . we see that g 
is a constant at a = 0, which then results in the fact that the first equation of (16. 6 p can be 



trivially integrated. But before doing that, let us find out which is the value of g\ 
that purpose, let us consider the normalization condition: 



1 



\3) 



With 



(6.9) 



where we have to integrate along (pi, 02 and an angle built in the "plane" of the two radial 
directions p, a (heuristically, think of introducing some polar coordinates r, 9 such that p = 
rsm9 and a = rcos9. Then we want to integrate in 9 from to | at fixed and large r). 
Inserting (16. 4p . we find 



g\{a=oo,p=0) — g\{a=0,p=oo) 

But from the first equation in (16. 6 p we read that g 



1 



p=0 



and, thus, g\ 



-1. 



now ready to integrate the first equation in (16. 6p at a = 0: 



3l 

47r2 



(^l<T=o) 



9ym 



47r2 L 



Nf 



Nf 



1\ f \ 1\ f 

Nc-^<S){p-Pq) log p+^e{p-pQ)\ogpQ + const 



where the next to last term comes from requiring continuity of the metric at p 
use of the radius-energy relation p = ^ found in |105j . we get: 



P{gYM){fj) 



gvM 

' 87r2 



Nf{p) ^ 



(6.10) 
We are 

, (6.11) 
Making 

(6.12) 



where Nf{p) is defined as the number of flavors for which the modulus of their masses is 
smaller than the scale. Matter fields with bigger mass are holomorphically decoupled at lower 
scales, as expected. The expression (I6.12p fits field theory expectations and is a non-trivial 
check of the described unquenched set-up. For further discussion of this model, see |127] . 



6.2 Flavors in lower-dimensional SQCD Models 

The approach described in subsection 16 . II can be also applied to construct supergravity duals 
of SQCD-like models in two and three dimensions by considering lower dimensional branes 
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wrapping different cycles of Calabi-Yau manifolds. In this subsection we will review two of 
such constructions. First of all, following refs. [661 1128] , we will consider the case of D3-branes 
wrapping a two-cycle of a Calabi-Yau twofold, which is dual to a two-dimensional gauge 
theory with M = (4, 4) supersymmetry. Secondly, we will review the similar construction 
of refs. [66| 11291 1112] of the gravity dual of three-dimensional N' = 4 gauge theories from 
D4-branes wrapping two-cycles in a CY2. Backgrounds dual to 2d and 3d flavored theories 
with reduced supersymmetry have been also constructed [521 113UI 1131] , and they will be also 
very briefly reviewed. 



6.2.1 Two-dimensional theories 



Let us consider the following setup for two sets of D3-branes in a Calabi-Yau cone of complex 
dimension two: 



CYa 









N2 




D3 (color) 
























Nf D3 (flavor) 























where S"^ represents the directions of a compact two-cycle and N2 are the directions of the 
corresponding normal bundle. Notice also that the symbols "— " and "■" represent respec- 
tively unwrapped worldvolume directions and transverse directions, while a circle denotes 
wrapped directions. Let us parameterize the cycle by means of two angular coordinates {9, (j)) 
and let a be the radial coordinate of the CY cone. The ansatz for the string frame metric 
which we will adopt is the following: 



ds 



St 



dxh + A; ( d9^ + sin^ 9 ( 



+ 



da'^ + a'^ \ dtp + cos 



+ dp^ + p^dVt 



where m is a constant with units of mass which, for convenience, we will take as: 



a 



(6.13) 



(6.14) 



Notice that in this setup there is another radial coordinate p, which represents the distance 
along W^, the directions orthogonal to both the D3-brane worldvolume and the CY cone. 
Moreover, is the metric of a unit three-sphere. Furthermore, the function z (which 
controls the size of the cycle) and the warp factor H should be considered as functions of 
the two radial variables (p, ex): H = H{p,a), z = z{p,a). 
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As in any background created by D3-branes, our solution should be endowed with a self- 
dual RR five- form F5, that we write as: 



F5 = + *J-5 . (6.15) 

The presence of Nj fiavor D3-branes induces a violation of the Bianchi identity of F5. Indeed, 
the WZ term of the fiavor brane action contains the term Y^Nf Im4 ^ ti^d.t acts as a source 
for this violation. Actually, the smearing procedure amounts to performing the following 
substitution in this term: 

J2 [ C'4 ^ / QACi, (6.16) 

where f2 is a six-form proportional to the volume form of the space transverse to the world- 
volume of the fiavor brane. The modified Bianchi identity takes the form dF^ = 2^1^ T3 Q. 
As in the four- dimensional example discussed in subsection 16.11 we shall locate the fiavor 
branes at a particular value p = pg of the p coordinate (the mass of the matter fields is just 
niQ = pQ/{2Tia')). Moreover, we will smear the Nf D3-branes along the angular directions 
{6, 0) of the cycle as well as along the external three-sphere. The corresponding smearing 
form is: 



^ = S{p - pq) dpAuj3Auj2 , (6.17) 



with 0J2 = sinOdO A d(f) and is the volume element of the external 5*^ with line element 
dQ^ (the minus sign in ( ]6.17p is due to the orientation of the worldvolume required by 
supersymmetry) . It is clear that the modified Bianchi identity in this case is: 

dF5 = -271 gs {a'Y Nf 6{p - pq) dp A UJ3 A U2 . (6. 18) 

Accordingly, let us represent F5 as in (16.151) with J^s being given by: 

J-5 = /s - 27ig,{a'fNfQ{p-pQ)cosAu2 , (6.19) 

with /s such that df^ = 0. We shall represent /s in terms of a potential C4 as /s = dC^, 
where C4 is given by the ansatz: 

C4 = gu^ A {dip + cos 9 d(j)) , g = g{p,a) . (6.20) 



Proceeding as in subsection 16. H one gets in this case the following set of BPS equations: 

i''[g-2^g,{a'fNfQ{p-pQ)] = p' z' , 



m 



m^H=—, g' = ~ap^H, 

a 

g = ^H' - ^Hm'[g - 2n g^a'f Nf Sip - pg)] , (6.21) 
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where the prime and the dot have the same meaning as in (16. 5p . The fulfillment of (I6.2ip 
ensures the preservation of eight supersymmetries by the background, which corresponds to 
M = (4,4) SUSY of the dual gauge theory. Moreover, one can prove that z{p,a) satisfies 
the following PDE: 



pziz - az) = a{pz^ + pz" + Sz' ) + ^ 5(p - p^) . (6.22) 

zJSc p"^ 



In the unfiavored case Nj = 0, the BPS system f l6.2ip (and the PDE equation fl6.22p ) can be 
solved analytically |128] by constructing the solution in five-dimensional gauged supergravity 
and by uplifting it to ten dimensions [66] . After a suitable change of variables one can show 
|128] that the metric and RR five-form of this gauged supergravity solution can be written as 
in our ansatz. In the general fiavored case one has to apply numerical techniques. However, 
as in the four- dimensional case, one only needs to know the solution for a = in order to get 
the behavior of the gauge coupling. Indeed, by means of a probe calculation one can check 
|128] that the supersymmetric locus of a color D3-brane occurs precisely at a = and that 
the gauge coupling is related to z{p, cr = 0) by means of the relation: 

^ -^(/^'^ = 0). (6.23) 



It follows from the system (16.211) that g{p,a = 0) is constant. Actually, by using a fiux 
quantization condition similar to the one employed for the 4d case, one can verify that 
g{p-,(J = 0) = where m is the constant defined in ( ]6.14p . By using this result in the 

first equation in (I6.2ip one readily integrates z{p,a = 0). By imposing continuity of the 
solution at p = pq, one gets: 



z{p,Q) = z, ^^NjQ{p-pq) iv, - -le(p-pQ) ,(6.24) 



P 



where z^, is a constant of integration. Plugging this result in (I6.24p . and assuming that the 
energy scale p is related to the holographic coordinate p as p = 27ra' p, one gets: 

9ym I .r Nf{p) 



(l-|^(-.-^)). (6.25) 



9YMip) 9ym ^ 27r/i 

where Nf{p) is again the number of fiavors with mass smaller that the scale p and qym is 
the bare UV Yang-Mills coupling. The dependence on the scale p of the Yang-Mills coupling 
displayed in (I6.25P matches precisely the one in field theory, which constitutes a non-trivial 
test of the gravity result. 

Backgrounds dual to 2d theories with M = (2, 2) SUSY can be obtained by wrapping 
D5-branes along a four-cycle of a Calabi-Yau threefold |130] . An alternative construction, 
which improves the UV behavior of the solution, involves D3-branes wrapping a two-cycle 
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of a CKj [661 1131] . One can further reduce the amount of supersymmetry by considering a 
D5-brane wrapping a four-cycle of a manifold of G2 holonomy, which leads to a dual of an 
M = (1, 1) supersymmetric gauge theory. In all these cases the flavor branes are extented 
along some of the non-compact normal directions of the cycle wrapped by the color branes 
and the corresponding backreacted solutions can be obtained numerically and are similar to 
the one reviewed here. 



6.2.2 Three-dimensional theories 



A similar analysis can be carried out to obtain the gravity dual of A/" = 4 three-dimensional 
gauge theories. In this case one must consider flavor and color D4-branes wrapping two-cycles 
according to the array: 



CY2 









N2 


M3 


D4 (color) 
























Nf D4 (flavor) 























The concrete ansatz for the ten-dimensional string frame metric we will adopt in this case 
is very similar to the 2d and 4d cases studied above, namely: 



dec ~\~ 



de^ + sin^ 



+ 



+e 



-2<i> 



\ida'^ + (di) + cos^rf^) j + dp^ + [dQ'^ + sin^ Qd^^) 



(6.26) 



where $ = $(p, cr) is the dilaton and the constant m is now given by: 

1 



SttQs Nc {a 



(6.27) 



As before, z = z{p, a) and the background should include an RR form whose Bianchi identity 
is violated due to the presence of flavor branes. For D4-branes the appropriate RR form is a 
four-form F4. If we locate the flavor branes at a fixed distance p = pg in the transverse 
and we smear them along their orthogonal angular directions, the modified Bianchi identity 
is: 



dF. 



2 ^]^Q -^"4 



2A', 8m3 



^{P - Pq) dp a L02 a L02 



(6.28) 



where U2 and 002 are the volume forms of the unit 
( I6.28P by means of the following ansatz: 



and 



two-spheres. Let us solve 



dCs + 



f 



2Nr 8m3 



0(P 



U}2 A U}2^ 



(6.29) 
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where C3 is the following potential depending on the function g{p, a): 

C3 = -guj2 A {d-^ + cosed^)). (6.30) 

By imposing that the system preserves eight supersymmetries, we arrive at the following 
system of BPS equations: 



Nf 1 



-4<I> 

e a 



1 



zz 



(6.31) 



2 -2 -4$ 

— m (7 2; e 



^7 + 



N 



f 



2Nr 8m3 



Again, one can combine the different equations in f l6.3ip and get a single second-order PDE 
for z{p,a), namely: 



z {z — a z) = pa {pz^ + p z" + 2z'^ + 



N 



a 



f 



1 



2Nr 8m3 



(6.32) 



As in the 4d and 2d cases, eqs. (16.311) and fl6.32p can be solved analytically when Nf = by 
using gauged supergravity [6611129] . In the general flavored case one can get analytically the 
form of the solution for a = [112] . Indeed, one can verify from fl6.3ip and the corresponding 
flux quantization condition that gip,^ = 0) = —l/{8m^) and that z{p,0) is: 



z 



(^'°) = ^*-8^2iV^®^^~^^^"^ 



1 - 



Nj_ 

2 Nr. 



0(p-pq) 



(6.33) 



with z^: being a constant. Moreover, by means of a probe calculation one readily verifles 
that 0" = is the SUSY locus of the color D4-branes and that the relation between the YM 
coupling and z{p,Q) is gyltip) — z{p,Q) / {2Tigs\f(y' m?). Using this result and the radius- 
energy relation p = 2na'p, one can convert (16.330 into the following equation for the running 
of the YM coupling of the 3d theories: 



1 



9YMiP) 



9ym 



1 - 



9ym 



Nr 



2 J 



(6.34) 



which again matches the corresponding fleld theory result. 

A gravity dual oi M = 2 three-dimensional gauge theory based on D5-branes wrapping 
a three-cycle was found in [1321 1133] . The addition of flavor to this background along 
the lines discussed here is carried out in [32l I131j . Alternatively, for this same amount of 
supersymmetry one can consider D4-branes wrapping a two-cycle of a Calabi-Yau threefold 
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7 A Mathematical Viewpoint 

In the approach we have followed up to now in this review on how to add unquenched 
flavor, we considered a family of equivalent embeddings of the flavor branes. This family can 
be generated by acting with the isometries of the background on a fiducial representative 
embedding. When the number Nf of flavor branes is large, considering the set of branes as a 
continuous distribution is a good approximation. We then computed the RR charge density 
generated by the branes, i.e. the smearing form Q, by explicitly performing the average over 
the set of embeddings and, subsequently, we have studied the deformation induced on the 
metric and forms due to the backreaction. 

The outcome of this microscopic approach is a system of supergravity plus delocalized 
sources preserving some amount of supersymmetry. It turns out that, in this process, very 
interesting mathematical structures emerge. The reason for this is the fact that the super- 
symmetric sources that we are using satisfy a calibration condition. As a consequence, one 
can use the methods of modern geometry to find backgrounds with smeared flavors in a 
systematic way [32j. In this procedure one does not deal with the set of embeddings and, for 
this reason, we will refer to it as the macroscopic approach, as opposed to the microscopic 
approach reviewed in previous sections. The goal is computing (or at least constraining) the 
smearing form Q by using the same type of technology as the one employed in the analysis 
of flux compactifications of string theory (see refs. |135] - [T38] . [29]). 

The central object in this geometric approach is the so-called 'calibration form' JC. For 
Dp-branes /C is a (p + l)-form, which can be represented in a vielbein basis as: 

r = i r (7 i) 

with e"^ ■■■"p+i = e*^! A ■ • ■ A e"''+^. The different components /Cai---ap+i are given by fermionic 
bilinears of the type: 

where e are Killing spinors of the background, conveniently normalized, and r is a constant 

matrix which (in the type IIB theory) is r = Tg ^ zr2, where T2 and T3 are Pauli matrices and 
the spinor e is represented as a two-dimensional vector of Majorana-Weyl spinors (r is the 
same matrix that appears in the expression of the kappa symmetry matrix of a Dp-brane 
when all worldvolume fluxes are switched off). The form /C can be used to characterize 
{p + l)-dimensional surfaces in the ten-dimensional geometry. A {p+ l)-dimensional surface 
Aip+i is said to be calibrated by K. if its pullback to Aip+i is equal to the induced volume 
form on A^p+i, namely: 

/C = v^^ditl rfP+^e , (7.3) 

where the ^'s are a set of local coordinates of Aip^i. When there are no NSNS fluxes or 
worldvolume gauge fields the calibration condition (17. 3p characterizes the supersymmetric 
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embeddings of Dp-branes (this can be easily established by using kappa symmetry). Actually, 
a Dp-brane whose worldvolume Aip+i is calibrated by /C is electrically charged with respect 
to an RR (p+2)-form field strength -Fp+2 and, in the Einstein frame, -Fp+2 is related to /C as: 

= d(^e"-^^IC^ . (7.4) 

Eq. (17. 4p is a consequence of supersymmetry |135j and, actually, in our backreacted back- 
grounds it follows from the system of BPS equations. Moreover, as a consequence of (17. 3p . 
the action of a localized embedding of a Dp-brane (without NSNS fiux and with worldvolume 
gauge fields switched off) can be written as: 



nloc rp 

•J Dp ~ -^P 



Mp+i 



(7.5) 



Following our prescription, the smeared version of the brane action is obtained by performing 
the wedge product with Vl of the (p + l)-form inside the brackets in (17. 5p and by integrating 
the result over the full ten-dimensional spacetime: 



nsmeared rp 

•J Dp ~ -'p 



e^^/C - Cp+^ 



A n . (7.6) 

'Mio '- -' 

Let us now define the (8 — p)-form -Fg-p, under which the Dp-brane is magnetically charged, 
as: 

Fs.p = ±6^^^ *Fp+2 , (7.7) 

where the sign depends on the particular value of p and on the conventions used. As in the 
examples studied in previous sections, the Dp-brane modifies the Bianchi identity of -Fs-p? 
namely: 

dFg.p = ±2nlQTpVL . (7.8) 

Eq. (17. 8p establishes a crucial connection between the smearing form f2 and the calibration 
form /C. Indeed, by using (17. 4p and (17. 7p . the right-hand side of (17. 8p can be written in 
terms of /C and its exterior derivative. Moreover, from the inspection of the smeared brane 
action (17. 6p . one concludes that Vl can be regarded as a kind of orthogonal complement (the 
Poincare dual) of K, in A^io- Interestingly, the possible calibration forms /C in a manifold 
are known and are related to its supersymmetric cycles and G-structures. In the case of a 
manifold preserving minimal SUSY in 4d, K, can be written in terms of powers of the Kahler 
form and of the holomorphic volume form. Thus, geometry and topology constrain the form 
of the charge density distribution of supersymmetric configurations and, actually, one could 
adopt the point of view in which the expression of VL is determined from these constraints 
without explicitly performing the average over the family of embeddings. This macroscopic 
approach was followed in refs. [32l 11341 1130t 1131] for some particular brane setups. 
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To finish this section let us detail the implementation of these mathematical concepts in 
the case discussed in section [21 namely the D3-D7 system. From now on we will assume that 
the metric, dilaton and forms are given by the expressions written in fl2.32p and f l2.33p . It is 
convenient to define the following two-form: 

J = [ e^f Jj,E + e^f dp A (rfr + Ake) ] , (7.9) 

which is such that h^^ J' is the Kahler form of the transverse 6d space. Actually, one can 
immediately verify that d[h~2 J''^ = as a consequence of the BPS equation for g in fl2.34p . 
By explicitly computing the fermion bilinear in (17. 2p and by using the projections satisfied 
by the Killing spinor of the flavored AdS^ x background, one gets that the calibration 
form /C in this case is given by: 

/C = ^ Vol(Mi,3) A J A J , (7.10) 

with Vol(Mi 3) = d'^x being the volume form of the Minkowski part of the space. Using 
the fact that dAxE = '2'Jke, one gets: 

d(e'^)C^ = ^ e^^^'^ (4^' + 0') e^s - 46^^ d^x A Jke A Jke Ndp = , (7.11) 

where, in the last step, we have used the condition (17. 4p for p = 7. Let us now verify that the 
value of Fg obtained in (17. lip is consistent with the expression for Fi written in our ansatz 
(I2.33P and, thus, with the displayed in eq. ( 12.3ip . Taking into account that the volume 
form of the KE space is \Jke /\ Jke, one can easily compute the Hodge dual of Fi and get 
the following result for Fg: 

Fg = -e^^ *Fi = ^ p{p) e^3+^^ d^x A Jke A Jke A dp . (7.12) 
The expressions (17. lip and (17.120 for Fg coincide if the following relation holds: 

Qfp{p) = e-^ [^g' + (p' - 4e2^-29] . (7.13) 
One can easily check that (I7.13P is a consequence of the BPS system (12.340 . 



8 Discussion 

In hindsight, we can say that the program of finding solutions dual to theories with un- 
quenched fundamentals with smeared flavor branes has been quite successful. As expected, 
it simplifies matters both when looking for the background solution and when discussing the 
physics they encode. 



90 



We have presented a series of example of solutions of ten-dimensional type IIA or type 
IIB supergravity coupled to D-brane sources. The philosophy and methods used in the dif- 
ferent cases are quite similar. Finding a consistent solution requires solving at the same 
time the closed string degrees of freedom (namely, finding solutions of the generalised Ein- 
stein equations in the presence of sources) and the open string degrees of freedom (namely, 
checking that the D-brane embeddings which generate the mass and charge source density 
are indeed solutions of the background). Supersymmetric solutions are easier to deal with 
and indeed, preserving SUSY simplifies enormously the technical work. It is rather remark- 
able that sometimes such complicated coupled systems can be (at least almost) completely 
integrated and the solutions given in a simple closed form (in particular in sections El |5l for 
the other sections, profuse numerical integration was necessary). However, supersymmetry 
is not mandatory for the construction and we have presented non-supersymmetric black hole 
solutions. 

The solutions are conjectured to be dual to theories with unquenched quarks. Since we 
have always dealt with the particular case of smearing the flavor branes over the transverse 
directions, we have built duals of a very particular class of such unquenched theories. We 
have used the solutions to discuss many physical features of the different set-ups. Many 
crosschecks of field theory expectations have been discussed. Just to mention a few instances, 
the running of the gauge coupling in different theories, the behaviour of the cascade in 
section [5] or the direct computation of the first flavor contribution to the entropy of the 
D3-D7 plasma (section 12. 5p which was previously known from an indirect method (namely, 
from first computing the free energy) [76]. All this asserts that the dualities discussed in 
this review are on firm ground. Due to obvious space constraints, we have not been able to 
include all the material that may deserve to be reviewed, but we hope that we have given 
enough references to the original literature. 

It is worth recapitulating about the presence of singularities in the different solutions 
discussed. First, in all the cases presented there are IR curvature singularities when all of 
the flavor branes reach the bottom of the geometry, see the heuristic picture of section 11.51 
They pertain to the kind of singularities usually called good. In fact, we have shown explicitly 
in the examples of sections |2] and [5] how adding (even an infinitesimal) quark mass leads to 
regular backgrounds (the analogous generalization for the set-ups of sections [3lll] remains an 
interesting open question). Moreover, heating up the theories can result in the formation of 
a black hole horizon behind which the IR singularity is hidden, see section [275] and eq. fl3.28p 
for examples. 

On the other hand, the solutions of sections [2] and [5] are singular in the UV (an effect 
connected to having flavor D7-branes) since the dilaton diverges at a finite position of the 
radial variable. This is expected, since it is the consequence of the Landau pole of the dual 
theory (more precisely, in the case of section [5] it is a duality wall). Despite the singularity, we 
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have shown that it is possible to consistently compute IR observables as long as the IR scale 
is well separated from the pathological UV. Clear examples are the meson spectrum (section 
I2.4p and the black hole thermodynamical properties (section 12.5. ip . In the D5D5 setups of 
sections [HI HJ the dilaton diverges linearly in the UV, signalling a little string theory-like UV 
completion of the dual field theory. We want to stress here that this already happens in the 
unfiavored solutions and thus is not problem associated to the backreaction. 

Finally, all the solutions in section [6] have a singularity in the IR. This singularity is not 
associated to the flavors as it is already present in the flavorless solutions and, at least in 
some cases, can be resolved by the worldsheet CFT [139] . On top of that, for the same 
models, typically, when Nf is sufficiently large, a Landau pole is generated and, jointly, a 
UV singularity appears in the geometry. 

We end this discussion with two clarifications. 

• We have repeatedly stressed that our main goal is to build duals to theories in which 
Nc and Nf are of the same order. Nevertheless, for the set-ups discussed in section [21 
which include the specially interesting flavored AdS^ x case, Nf -C Nc is needed, 
see section [2^6] (similar comments apply to section [5]). This is because, starting with a 
conformal theory, the introduction of extra matter generates an UV pathology, namely 
a Landau pole. Then, roughly speaking, in order to have a meaningful IR, it has to 
be well separated from the pathological region, enforcing the number of flavors not to 
be too large. However, backreaction effects and, accordingly, the effect of unquenched 
quarks, can still be computed as an expansion in Nf/Nc. On the other hand, for the 
models in sections [3l [H [6l this restriction is not present and, indeed, it makes sense 
to talk about solutions with Nf ~ Nc- In fact, this is imperative for instance when 
discussing Seiberg-like dualities as in section 13.6.21 

• Since the (DBI) action is used to model the D-brane sources, one could be wary for 
the following reason: the effective string coupling on a stack of Nf D-branes is QgNf ~ 
Nf/Nc and this should be small for the DBI to be a good approximation |140] . whereas 
there is not a good effective description for strong string coupling. However, this caveat 
is circumvented because we do not deal with stacks of localized flavor branes: due to 
the smearing, the typical distance between any pair of flavor branes is of the order 
of the size of the transverse space, which is typically large. As a result, the flavor 
symmetry is usually broken to f/(l)^^ and the effective open string coupling remains 
small. As already pointed out, this amounts to keeping "one window graphs" in the 
Veneziano expansion [25]. 
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Outlook 



There are still many open questions that deserve to be addressed within the framework 
presented in this review. We briefly mention a few examples of possible future projects. 
They comprise both making further progress in studying the models here presented and 
building new solutions that could be useful in exploring the consequences of the formalism 
for different physical points. Along the first of these lines, it would be nice to generalise the 
solutions of section |3] to the massive quark case in order to remove the IR singularity. Also, 
we expect the black hole solution of section 12.51 to encode interesting physical information. 
For instance, one could consider massive embeddings in the search of a first order phase 
transition similar to those in [T6] . 1141] . [76]. The peculiarity of the back- reacted setting 
would be, conceivably, that the area of the horizon would undergo a finite jump at the 
transition. Along the second line, a back-reacted D4-D6 solution building on the model of 
[16] could be useful in discussing QCD-like properties. Another conceivable program is to 
look for a solution, which, as in [142] . may correspond to a color- flavor locking phase. The 
study of fluctuations in these backgrounds, that will also contain fluctuations of the fields 
in the flavor branes, with a view on understanding holographic renormalization would be a 
highly interesting result. 

Aside from this, it would be nice to find solutions (with backreacted flavor branes) that 
contain an AdS^ factor. The study of conformal anomalies there may give interesting results. 

As stressed in the introduction, finding the kind of solutions discussed here, including the 
D-brane backreaction, has an interest on their own, independently of AdS/CFT. It would be 
nice to understand whether they may turn out to be useful for different physical applications. 
For instance, for models of inflation built with D3-D7 systems on the conifold (see [143] for 
recent progress in this direction), the analysis of section |5] could have some relevance. 
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